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Abstract 

We study linear time dispersive and dissipative systems. Very 
often such systems are not conservative and the standard spectral 
theory can not be applied. We develop a mathematically consis- 
tent framework allowing (i) to constructively determine if a given 
time dispersive system can be extended to a conservative one; (ii) to 
construct that very conservative system - which we show is essen- 
tially unique. We illustrate the method by applying it to the spectral 
analysis of time dispersive dielectrics and the damped oscillator with 
retarded friction. In particular, we obtain a conservative extension 
of the Maxwell equations which is equivalent to the original Maxwell 
equations for a dispersive and lossy dielectric medium. 
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1 Introduction 

In this paper we describe a mathematical framework for a spectral the- 
ory of linear time dispersive and dissipativc (lossy) media, e.g., dielectric 
media. Here is a concise formulation of the setup. We consider a linear 
system (medium) whose state is described by a time dependent generalized 
velocity v (t) taking values in a Hilbert space Hq with scalar product (•, •). 
The evolution of v is governed by a linear equation incorporating retarded 
friction 

mdtv (t) = -lAv (t) - / a (t) v {t - t) dr + f (t) , (1.1) 

where to > is a positive mass operator in Hq, A is a self-adjoint operator 
in Hq, f (t) is a time dependent external generalized force, and a (t), t > 0, 
is an operator valued function which we call the operator valued friction 
retardation function |H| Section 1.6], or just friction function. The names 
"generalized velocity" and "generalized force" are justified when we inter- 
pret the real part of the scalar product Re {(w (t) , f {t})} as the work done 
by / (t) per unit of time at instant t. Hence, the total work W done by 
fit) is 

W= / Re{{v{t),f{t))} dt. (1.2) 
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By the same token, the first two terms on the right hand side of 
namely 

/•oo 

-iAv{t)- a{T)v{t-T) dr, (1.3) 
Jo 

are interpreted as the force which v exerts on itself at time t. This self forc- 
ing depends on v through two terms: (i) the instantaneous term ~iAv {t) 
and (ii) the "time dispersive" integral term a (r) v {t ~ t) dr involving 
values V {t') for t' < t. The integral term is interpreted as a (retarded) 
friction force whose special form reflects two fundamental requirements: (i) 
time homogeneity and (ii) causality. 

If we rescale the variables according to the formulas 

V = ^/mv, il = ^=A^=, a = ^=a^=, / = ^=/, (1.4) 



then the equation Hl.l|) reduces to the special form with m the identity 
operator, i.e. 

/"OO 

dti (t) = (t) - h{T)v{t-T) dT + f{t). (1.5) 

We refer to Q. — -^=A^^ as the frequency operator in Hq since its spec- 
trum gives the resonant frequencies of the system. 

We call the system conservative, or non-dispersive, if the friction func- 
tion vanishes, i.e., 

mdtv (t) = -lAv (t) + fit). (1.6) 

In this case, it is natural to define the internal energy of the system at time 
t to be ^ {v{t),mv{t)) , since it follows from (|1.1() that 

~ivit),mvit))^Re{{vit)jm, (1-7) 

where Re {{v (t) , f {t))} is the instantaneous rate of work (by assump- 
tion). In the non-conservative (dispersive) case, we continue to interpret 
^ {v{t),mv{t)) as the internal energy but note that H1.7|l is not true in 
general. Instead, 

~ {v{t),mv{t)) 

^Re{iv{t)J{t))}-Rc{(v{t)J aiT)v{t-T) dr]} . (1.8) 
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We interpret the second term of this expression as the instantaneous rate 
of "work done by the system on itself," or more properly the negative rate 
of energy dissipation due to friction. 

In all physical models of which we are aware the time dispersive term a 
arises as a phenomenological description of a linear coupling between the 
system and some other "hidden" degrees of freedom. Our belief in the 
conservation of energy suggests an equivalent description with the hidden 
system described by a vector w in a (different) Hilbert space Hi, such that 
the extension 



evolves according to a conservative equation like H1.6|l . 

The main result of this paper is that under a physically reasonable 
hypothesis it is possible to construct a conservative extension of the original 
system (II. 1|) in the following form 

mdtv (t) = -iAv (t) - iTw (t) + f {t) , to > 0, A is self-adjoint, (1.10) 
niidtw (t) = —iT^v (t) — lAiw (t) , toi > 0, Ai is self-adjoint, (l-H) 

where w ^ Hi, the Hilbert space of "hidden" variables, mi and Ai describe 
respectively the mass operator and the generator of internal dynamics on 
Hi, and T : Hi Hq is a coupling operator between the "hidden" and 
"observable" variables. While it may be conceptually clear that such ex- 
tensions should exist for reasonable linear models, it is not immediately 
obvious how to determine from the given equation Hl.lj) if such an exten- 
sion does exist. Nonetheless, we show there is a natural intrinsic condition 
- non-positivity of the total work done by the friction force - which is both 
necessary and sufficient for such an extension. Furthermore, we shall see 
that the Hilbert space Hi of "hidden variables" as well as the operators 
Ui = -^Ai-j^ and r in H1.10|) - (|l.ll|l are essentially uniquely determined 




by the friction function. 

It would be interesting and natural to study (|1.1U|) . (|1.11|) . and thus 
p.l|l . with a random force f(t) as a fluctuation-dissipation model similar 
to the Langevin equation 8 . However, the importance of the subject 
and the efforts needed to conduct such a study are worthy of a separate 
publication. 

For a homogeneous dielectric medium described by a scalar-valued fre- 
quency dependent electric susceptibility x(w), and magnetic permeability 
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11 = 1, the conservative extension takes the form (in common notations) 

9tH(r,t)--VxE(r,i), 
5tE (r, t) = V X H (r, t) 



f ^JruxUx (<T)*(r, i,cr) da - AttJ {r,t) , (l-^^) 



m^dt'^ (r, t, a) = —im^a'^ (r, t,a) + m<^n<^ (f )E (r, t) , 

where the field ^ (r, t, a) describes "hidden" variables, which one may view 
as a "string of dipoles" with string coordinate a G R attached at every 
space point r. The parameters n-^ (a) and rriy^ are related to the electric 
susceptibility x{u!) as follows 

/•oo 

(f^) = 41m{CTx (a)} > 0, mT^ = / (cr) da. 



.CO . / (1-13) 

47rwx(^^)= lim / rfo", ImC > 0. 

V^+'^J-oo a - UJ - IT] 



In the second equation of H1.12() . the term yj'm<^n<^ (o')'S' (r, i, cr) dcr is 
47rx the displacement current, and thus the electric polarization is related 
to * by 

P(r,t) = -^/" I ^vfi^n^ (a)*(r,r,a) dadr, (1.14) 

assuming P = at t = — oo. In addition, the fields H(r, i), D(r,t) = 
E (r, <) + P (r, t) are required to be divergence free, i.e. 

V-H(r,t) = 0, 

V-E(r,t)+/ / ^msAi^ (cr)V • * (r, r, a) dadr = 0. (1.15) 

A system of extended Maxwell equations similar to (|1.12l) were proposed 
in [TH|. 

In many physical models there is an instantaneous contribution to the 
friction. Thus we assume throughout that a (t) , f > is of the form 

a (t) = aooS (t) + a{t) , where a (t) is continuous for t > 0, (1-16) 

with Qfoo self-adjoint (any anti self-adjoint piece can be incorporated in 
—iA). The representation H1.16|l explicitly distinguishes the instantaneous 
component aooS {t) of the friction fmiction from the retarded component 
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described by a (£) . It is possible to extend some of the results below to 
a(t) which are operator valued measures or distributions, but to simplify 
the exposition such extensions - which are physically somewhat esoteric 
anyway - will not be considered. On the other hand, the instantaneous 
component aoo^ {t) is completely natural and occurs in many examples. 

To state the central condition of this paper, let us consider the total 
work Wfr done by the friction force, which may be written as follows 



( a{t) if < > 
ae{t)^2a^6{t) + l Rc{a(+0)} if < = , - oo < i < cx). (1.18) 
[ at {-t) if < < 



Notice that ae{-t) = al(t). The precise condition which distinguishes 
systems with conservative extensions is that Wfr should he non-positive, 
that is. 



for every function v{t), signifying the ultimate conversion of mechanical 
energy into heat ~ the transport of energy from relevant (observable) de- 
grees of freedom to hidden ones. We refer to (|1.19|l as the power dissipation 
condition. 

Scalar functions (t) which satisfy (|1.19() . known as positive definite 
functions, are familiar from Bochner's theorem on the Fourier transform of 
a finite positive measure. In fact, one construction of a conservative ex- 
tension to (|1.1|) is based on Theorem l3.2l an operator valued generalization 
of Bochner's theorem. 

It is often useful to work in the frequency domain, and the power dis- 
sipation condition 11.19|l may be formulated there as well. To do so we 
formally apply the Fourier transform defined by 




(1.17) 



where 




(1.19) 




(1.20) 



to recast the time evolution equation (|1.1|) as 



wmv (lo) = [A — ia (lo)] v (lu) + if (lo) , 



(1.21) 



Here 




(1.22) 
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is a formal object, and may not be defined pointwise since we do not assume 
integrability of air). However, if it is defined pointwise, we see - in view of 
Bochner's theorem - that the power dissipation condition H1.19|l becomes 

Oe (w) > for all w, (1.23) 

where 

Oe (w) = 2 Re {a (w)} a (w) + (uj) . (1.24) 
That is, (Xg (lu) is a positive semi-definite operator, 

(w,ae (w)w) > for all u e 7i. (1.25) 

We refer to (|1.23|) also as the power dissipation condition. If the Fourier 
transform is not defined pointwise the situation is a bit more delicate, but 
in effect (to) duj is a non-negative operator valued measure (see Theorem 

roi l 

In many examples a{t) is Hermitian for every t. In those cases the real 
and imaginary parts of d (w). 

Re {d {lu)} = ^ [a (w) + (w)] , Im {d {cu)} = ^ [« (w) - (w)] , (1.26) 

are respectively even and odd functions of lu given by the cos and sin 
transforms of a. 



He {d {lu)} = a oc + / a {t) cos {lut) dr, 
Im{a(cj)}= / a (r) sin (cjt) dr. 



In the important special case of a conservative system p.6|l . equation 
H1.21|l reduces to 

Lumv {uj) = Av {i^) + if {u) . (1.28) 

When the external force vanishes (/ = 0), a time harmonic solution v {t) = 
V^e~'"* with V {uj') = V^5{lu — lu') is obtained for which solves the 
spectral problem 

LumV^=AV^. (1.29) 

This eigenvalue problem is part of the standard spectral theory of self- 
adjoint operators, and its analysis is instrumental to the study of the non- 
dispersive evolution (|1.6|l . In particular, we remind the reader that the 
eigenvectors VL, - which may lie in a proper extension of the Hilbert space 
- are m-orthogonal for different lo and form a basis of the Hilbert space Ti.. 
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In contrast, a time harmonic solution v (t) = Wi^e""^* to satisfies 

umv^ = A (lu) Vi^, with A {lj) — A — la (w) . (1.30) 

In non-trivial examples, the operator A [uj) is typically non self-adjoint. 
Consequently the spectral theory for A [lj) may be rather complicated, 
even in the simplest examples, with a finite dimensional Hilbert space, m 
the identity matrix, and A (w) a finite square matrix. For instance, A [uj) 
may not be diagonalizable as it may have nontrivial blocks in Jordan form. 
Consequently, the genuine eigenvectors of A [lj) may not form a basis. 
In many problems of interest, including continuum dielectric media, the 
Hilbert space is infinite dimensional and the operator A {uj) is unbounded, 
in addition to being non self-adjoint. Thus, in general it seems to be very 
difficult to analyze the eigenvalue problem H1.30|l . There are some results 
on the completeness of so called root vectors for dissipative linear opera- 
tors, Chapter V], but the conditions of those statements - for instance, 
compactness of the imaginary part of the operator - are too restrictive and 
are often not satisfied in problems of interest. 

Compounding these difficulties, even if the complicated spectral analy- 
sis of non-self-adjoint A (w) were somehow addressed for fixed w, another 
problem would arise since A {lj) depends on the spectral parameter uj. Thus 
there is no obvious relation between solutions to H1.30|l for different values 
of UJ. In particular, the vectors v^^ need not be orthogonal and it is not 
clear that they span the Hilbert space. 

However, for physically meaningful examples A {uj) is not an arbitrary uj 
dependent non self-adjoint operator, but rather one with certain properties 
which allow some kind of spectral analysis. The relevant properties, as we 
will see below, are analyticity and a dissipation condition for complex oj in 
the upper half-plane. 

The focus of this article is the construction of a spectral theory for 
a wide class of dispersive and dissipative systems - described by Hl.l|l . 
(|1.21|l . or H2.37|l below - under the assumption of a suitable power dissi- 
pation condition. As indicated above, our approach is based on the obser- 
vation that dispersion and dissipation are caused in most physical models 
by coupling with degrees of freedom which are "hidden" from observation. 
Consequently, we consider a spectral theory for a dispersive/dissipative 
system to be a realization of the system as a proper projection of a con- 
servative extension to which the standard spectral theory may be applied. 
As we show, the power dissipation condition is necessary and sufficient for 
such a conservative extension and determines the (minimal) conservative 
extension uniquely. 

The paper is organized as follows. In section 2 we recall the analysis of 
a general conservative system of the form (II.IOII . Hl.ll(l . which we take as 
an abstract model for a system with "hidden variables." In particular we 
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show that positive power dissipation holds for any system with a conser- 
vative extension. We also describe the admittance operator reformulation 
of and (|1.1(J|) . In section 3, we discuss the main mathematical 

results which demonstrate that positive power dissipation is equivalent to 
the existence of a conservative extension. There are two approaches to con- 
structing a conservative extension, proceeding by either the time or 
the frequency H1.21|l representations and based on the classical Bochner's 
Theorem 13.11 and Herglotz-Nevanlinna Theorems 13.81 13.91 respectively. In 
Section 4, we summarize in concise form the schemes by which one may 
construct a conservative extention of a given dispersive/dissipative system 
including (i) the space of "hidden" variables; (ii) a self-adjoint operator 
describing the internal dynamics of the "hidden" variables; (iii) an oper- 
ator coupling the "hidden" variables with the "observable" variables. In 
Section 5, we give examples of conservative extensions for general scalar 
dispersive systems and homogeneous isotropic dielectrics, obtaining an ex- 
tended conservative system which is equivalent to Maxwell's equations for 
a dispersive and lossy dielectric medium. In forthcoming work we shall de- 
scribe in greater detail the application of the techniques presented here to 
dispersion in dielectric media. In section 6, we discuss dissipation, or loss 
of energy, a phenomenon which can arise in systems with dispersion. We 
present sufficient conditions on a {t) for solutions v (t) to to exhibit 

dissipation, that is limj^oo 1^(^)11 = 0. Finally, Section 7 is devoted to 
proofs of results in Section 3 and related constructions. 

2 Modelling hidden degrees of freedom 

In this section we introduce and analyze abstract models for hidden degrees 
of freedom, showing in particular that the power dissipation condition holds 
for the truncation of any conservative system. In addition, we recall the 
equivalent description of a linear system in terms of its admittance operator 
and discuss the admittance for truncations. 

2.1 Abstract model: preliminary analysis 

Consider a conservative system with degrees of freedom (variables) divided 
into two classes: the observable variables, denoted v, and the hidden vari- 
ables, denoted w. For instance, v might describe the polarization of a 
dielectric medium at different space points, while components of w account 
for microscopic degrees of freedom which give rise to the material relations. 
We assume that v and w take values respectively in Hilbert spaces Hq and 
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Hi, and that the combination 



V 



V 

w 



eH = Ho® Hi, V e Ho, w e Hi, 



(2.1) 



describes a conservative system governed by the following evolution equation 

MdtV (t) = -lAV (t) + F (t) , (2.2) 

where A4 > and A are self-adjoint operators in Ti.. Furthermore we 
assume - and this is key - that the external force F (t) is of the form 



Fit) 



fit) 




= :Tt/(t) 



(2.3) 



indicating that we may excite directly only the degrees of freedom corre- 
sponding to the observable variables v. 

In H2.3|l we have introduced the operator: 



(2.4) 



where Ihq a-nd are respectively the identity operator in Ho and the 
orthogonal projection onto ran 

= Fo © {0} in n . Notice that : 
Ho — *■ 7i is an isometric injection from Ho into Ti., and we can recast the 
equations (|2.2|l . (|2.3|) as 



MdtV (t) = ~iAV it) + T^f (<) ,V{t)eH,f (<) e Ho. 
The operator T is an example of an isometric truncation: 



(2.5) 



Definition 2.1 Given two Hilbert spaces H and Ho, a bounded linear 
operator T : Ti Ho is called an isometric truncation of Ti. to Ho if 

Observe that when Hq C H the orthogonal projection Phq of H onto Hq 
is an isometric truncation. The conservative system as described by the 
equation ^2. 5p with M > 0, A self adjoint, and T an isometric truncation 
is of the most general form we consider in this paper. 

We are specifically interested in conservative extensions to (|l.l|l of the 
form (|1.10() - (|1.11() . corresponding to the special case of equation (|2.5I) with 
operators M and A in the following block-matrix form 



M = 



m 






A^ 



, m > 0, 

A r 
rt ill 



= A, n\^ni, r : Hi ^ Ho, (2.6) 



10 



where Ihi is the identity operator in Hi. In physical models with un- 
bounded operators, the block decomposition H2.6|l may be somewhat for- 
mal, since there remains the question of specifying the domain of each 
operator. In general we consider A and f^i which are self adjoint with 
domains and V^rii) respectively. For the time being, we take the 

operator T - which evidently provides the coupling between the observable 
and hidden degrees of freedom - to be a bounded map from Hi to Hq. 
In this case, the natural domain for A is the direct sum "D^A) © T>{fti), 
and with this choice the operator is self-adjoint. We consider unbounded F 
below to describe systems with non- vanishing instantaneous friction, e.g., 
a damped oscillator. In that case, F will be a map from T) (Sli) to Hq, and 
we shall have to specify the domain of A carefully. 

We can represent a formal solution to 12.2(1 - using the rescaling (|1.4() 
- by the formula: 

/•oo 

Jo 

= M-h''-^'"^'-'''^M-^F{T) dr, Am = M-^AM~K (2.7) 

assuming the system was at rest with V{t) — 0, F{t) = in distant past. 
In particular, we note that 



V {t) = M^h-'-^^^^M^^Fo, t > 0, 



(2.8) 



for a pulse force F (t) — FqS (t), and 



V{t)=M 



M~^-F{t) rfT + 0(l) 



(2.9) 



for i ^ oo if, say, M. ^F{t) < oo. 

To obtain an effective equation for the evolution of v, let us look at the 
block form of the evolution equation ((2.2|l with M., A and F respectively 
satisfying H2.6|l and (|2.3II . 



mdtv (t) = -iAv {t) - iTw (t) + / (i) , 
dtw (t) = -iF^w (t) - iniw (t) . 



(2.10) 
(2.11) 



Note that — iF^^u (t) plays the role of an external force in eq. H2.11|l . Thus 
we can solve for w as in eq. (|2.7|l . obtaining 



w{t) 



-1 / e-'^^'^VK {t - t) dr. 
'o 



(2.12) 
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Plugging this into cq. I|2.1U|1 yields 

mdtv {t) ^ -iAv {t) - Te-'^^^^T^v{t^T) dT + f{t). (2.13) 

"'0 

The dispersive evolution equation 12.13|l is of the form ll.l|l with friction 
function 

a(t) = re"'"i*r^ for t > 0, (2.14) 

and describes the evolution of the observable variable v. We note in par- 
ticular that the instantaneous friction Uao vanishes - since a(0) TT^ IS 
finite - and the extended friction function Og (t) is 

fle (t) = re-'^^i*r^ for -oo<t<OQ. (2.15) 

Notice that for ae (t) of this form we have 



OO pQO 



{v (t) ,a^{t — t)v (t)) dtdr 



OQ J —OO 



jifiltpt^ (i) dt 



(2.16) 



readily implying that the extended friction function (t) — Fe'^^^'^r^^ sat- 
isfies the power dissipation condition H1.19|l . The central point of this work 
is that f2i and T needed to satisfy i2.14{ l can be re- constructed from the 
friction function a{t). 

2.2 Unbounded coupling and instantaneous friction 

Let us briefly consider how to define the abstract model H2.6|l with un- 
bounded r, a necessary step in the description of systems with non- vanishing 
instantaneous friction (aoo 7^ in (|1.16|) ). The following discussion is by 
nature a bit technical; we direct the reader to the case of the damped oscil- 
lator in Section 5 for an explicit example which may provide clarification. 
Suppose r is given as a map T : I?(ili) Hq which is ili-bounded, i.e. 

Wm' <c\\{ni + ii„,)^\ (2.17) 

where Ihi is the identity operator and C < 00. We continue to take A 
and fii to be self adjoint on their respective domains ^{A) and X'(Oi) as 
above. The operator r<i>i^, with 

$fi = , for i? > 0, (2.18) 

is a bounded map from Hi to Hq, with ||r$fl|| < i? as i? ^ 00. We denote 
the adjoint of this map by $_Rr^, although has not been defined,^ and 

^It may be defined as a map from Hq to a proper extension of Hi, namely the space 
"D (B)* of conjugate Unear functionals on "D (B), but we do not use this fact here. 
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define A as the limit 



A 


V 


:— lim 




w 


_R — ^oo 



Av + T^rw 



(2.19) 



on the domain T> (A) of vectors [v, w]^ such that the hmit exits. 

As things stand, it is not clear if the resulting operator is self-adjoint, or 
even that V (A) is dense. To proceed we require an additional assumption 
- H2.24|l below - which guarantees self-adjointness. To state that condition 
consider the map 

s -.^Tini + iiH.r', 

which is bounded and satisfies 



(fli - UhJ^rS^v = <i>flr^w for any v G Hq. 
Therefore the limit 



(2.20) 
(2.21) 



lim ^Rr\' + rii^RW 



exists if and only if 



lim (ill 



+ S''v] (2.22) 



(2.23) 



which is equivalent to saying that w = cf) — S^v with </> G 'D{ili). Thus 



[v,w] G I? (-4) if and only if w = 

Gv lim r$i^S'^w = hm T 



S^v with G X'(fli) and 
1 



(2.24) 



exists. We require of T and fli that the limit H2.24|l exists for every v G 
V {A) , and defines an A bounded operator G with A bound less than one 
- i.e., there are i5 < 1 and Ps > such that 



\Gv\\ <6\\iA + iPslHo)v\\ 



Under these assumptions 
V(A)-- 



V 

w 



veV{A) and lim ^flT'l'v + fli<I>_RW exists i 
V eV{A) and G V {Vli) 
and the operator A is self-adjoint: 



V 

4> — Sv 



(2.25) 



(2.26) 
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Proposition 2.2 Let A, fii be self adjoint operators on the Hilbert spaces 
Hq, Hi with domains 2? (A), T) (f^i) respectively. Suppose that F : V (fii) 
Hq is such that (|TT7|l holds and the limit (|^!^ exists for aU v € V{A), 
defining an operator G for which the bound (|2.25|) holds. If A is defined 
by lfTTO|l on the domain V (A) C Hq Q) Hi specified in then A is 

self-adjoint. 

The proof of this proposition is elementary. It is obvious that A is 
symmetric, so to prove self-adjointness we need only to show that A^V ^ 
±iV implies ^ = 0, which is an easy exercise. 

For the operator A defined in this way, the evolution equations H2.1(J|I - 
H2.11|l imply that 

dt^RW (t) = -i<^Rr^v (t) - ini<^>RW {t) , (2.27) 

and thus 

{t) = -i / e-'*^^i^$flF^w {t - t) dr. (2.28) 

Therefore 

/•oo 

mdtv [t) ^ -iAv [t) - lim / Fe-'^i^^^F^u (t - r) dr + / (t) . (2.29) 

Jo 

While this equation is formally similar to H2.13|l . we note that the re- 
sulting friction function 

g-iOit 

ae{t) = lim F--3 ^^ for - oo < i < cx3 (2.30) 

Rl ^ ^Hi 

is defined only as a distribution, 

-init 

/ F^ 

R — ^oo 



ae{t)v{t) dt ^ lim / F-^^^ T'<v{t)dt, 



for veCc (M, Ho) , (2.31) 

and may not in fact be a function. In particular, there may be non- 
vanishing instantaneous friction. However a^it) is a relatively tame distri- 
bution; it may be expressed as a second order differential operator applied 
to a (strongly) continuous function 

ae{t)= + Se-'^^^'S'' with S = T{ni+ilH,y\ (2.32) 
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In a key example, Hi = L'^ Hq), the space of square integrable 
-ffo-valued functions on the real line, and Qi is multiplication by the inde- 
pendent variable, fll^p (x) = xip (x). Given any positive operator aoo, say 
bounded, on Hq we define 

[T^] {x) = —= / ^/a^ip (x) dx, (2.33) 

which is f2i bounded since 

1 f°° 1 

[Sij] (x)^^ —^V^i^ (x) dx (2.34) 

is bounded. Thus 

^g-ifiit^t ^ _ / dx^ aoo = e-l*laoo, (2.35) 

TT J-oo ^ + 1 

and ^ 

ae(i) = (^-^ + 5e-'"i*5^ = 2ac,5 (i) , (2.36) 

purely instantaneous friction. We shall return to this example in our dis- 
cussion of the damped oscillator below. 

2.3 Linear response and the admittance operator 

In the linear response theory, e.g. ref. Section 3], a system is often 
characterized by its admittance operator 21 (oj) : Hq Hq, defined by the 
relation 

v{uj)^^ (lu) f (uj) , v{lu), f (uj) e Ho (2.37) 

at real frequencies uj. Here we discuss the admittance formulation of dissi- 
pative dispersive systems satisfying the power dissipation condition (|1.19|) 
and truncated conservative systems described by H2.5|l . For this purpose, 
it is useful to recast the linear response equation H2.37I) in the domain of 
complex frequencies C = ^ + i??, ?? > 0, which corresponds to replacing the 
Fourier transform by a Fourier-Laplace transform. 

We assume that the system governed by (11.1(1 is at rest for all negative 
times, i.e. 

v{t) = 0, f{t)^0, t<0, (2.38) 
and define the Laplace transforms 

HO^ e'^'v{t)dt, f (0 = e<'f{t)dt, (2.39) 
Jo Jo 

a (C) = / e'^*a {t) dt ^ + [ e"^*a {t) dt, (2.40) 

J[0,oo) J[0,oc) 
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on the domain {Im^ > 0}. It is convenient to assume that a {t) is bounded 
as t ^ oo so that a (C) is well defined (see Condition 13.51 below). Under 
the Laplace transform, the evolution equation 1)1. l|l is transformed into the 
following identity 

Cm« (C) = [A-m (C)] « (C) + i/ (C) , C = + i??, ^ = Im C> 0. (2.41) 
We have Re a (C) > (for Im^ > 0), since 

Re (w, a (C) v) 

^ ^ (w, Oe (t - r) w) e'«*e-'f^dtdT > (2.42) 



[0,oo) J[0,oo) 

for any v E Hq by the power dissipation condition (|1.19() . Hence, the 
operator — yl + ia {() is invertible, and 

V (C) = ^m.A.a (C) / (C) , C = ^ + = Ini C > 0, (2.43) 
2l™,A,a (C) = i[Cm-A + ia (C)]"' . (2.44) 

The equation H2.43|l generalizes (|2.37|) to a certain extent since it is an 
identity for analytic functions in the upper half plane. Note that 

Re {^mA.a (C)} = I ["^mA.a (C) + <^A,a (O] 45) 

= 2l™,A,a(C){ImCm + Rea(C)}<,^,,(C) > 0. 

which expresses the power dissipation condition p.l9|l in terms of the ad- 
mittance operator 21^, a, a- 

The admittance equation (I2.43|) provides an essentially equivalent de- 
scription of the system Hl.l|) . In particular, the various operators in (|l.l|l 
can be readily recovered from 2lm.A,a (^) by the relations 

m"^ = - lim r]%ri,A,a i^v) , A = - lim ^/mlm'^lZ^A „ (iv) V^, 

a (C) = i {Cm -A) + [%n,A,a iOr' ■ (2.46) 

Hence, the admittance operator 21m, a. a i^) carries all the information about 
the system initially described by the triplet {m, A,a (()}. Very often the 
admittance equation H2.37|l is a preferred form, since the admittance 21 (w) 
may be measured experimentally more readily than to, A or a{uj). 

There are also several technical advantages to the admittance formula- 
tion l|2.43|l . First, the quantities v d {(), 21 {(), f {() are analytic func- 
tions in the upper half-plane Im^ > 0, whereas their time counterparts 
may be more singular functions. In addition (|2.43() has the advantage that 
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auxiliary operators which appear in the analysis are bounded. In particu- 
lar, we shall see that for aoo 7^ the admittance formulation permits us to 
avoid the subtleties required for an unbounded coupling F. 

Consider now a system which is the truncation of a general conservative 
system of the form (|2.5|l . Under the Laplace transform (|2.5|) is transformed 
into 

-iCMV (C) = -iAV (C) + (C) , Im C > (2.47) 
which is easily solved for V, 

ViO^HCM-Ay'T^fiC), ImC>0. (2.48) 
by the isometric truncation T :'H ^ 



Multiplying of the both sides of 
Hi) yields 

«(C) =21(0/(0, C = ^ + "?, r7 = ImC>0, (2.49) 

%{C)^iT{CM-Ay^T\ (2.50) 

Notice also that the admittance operator as defined in (|2.49l) satisfies 



Rc{2l(0} = 



21(0 + at (0 



it 



(2.51) 



^{lxnC)T{C,M-A)~^ M T{C,M-Ay^ 
This identity together with M > implies 

Re{2l(0}>0. (2.52) 



For operators M.,A\n the block- matrix form (|2.6|) we have 

2t (0 = i (c™ - ^ - r {CIh, - ni)-^ Ft) , (2.53) 

as may be easily verified. Comparing (|2.53|) and H2.44|) gives the follow- 
ing formula for the Laplace transform a (0 of the friction function for a 
truncated conservative system 

a{c) = ir{ciH,-n,y'r\ (2.54) 

which could also be verified by directly transforming (|2.14|l . 

A second approach to the central construction of this work is based 
upon the observation that the Hilbert space Ti and the operator triple 
Ai, A, T in eq. (|2.50|) can be reconstructed from the admittance 21 (0 • 
Alternatively -ffi, F, and fii appearing in H2.54|l can be reconstructed from 
aiC). 



17 



3 Bases for a conservative extension 



We now describe how, given a dispersive system in the form or H1.21|l . 
one can reconstruct the hidden degrees of freedom. Of course, the resulting 
mathematical reconstruction is initially devoid of physical interpretation. 
However, the reader should bear in mind that one usually knows a given 
evolution equation of the form involves "hidden" degrees of freedom 
with a natural physical interpretation. Generally, the abstract extension 
may be interpreted therefore in a physically concrete way. 

For example in a classical dielectric medium, the time dispersion comes 
from the material relation between the electric displacement D (r, t) at a 
point r and the electric field E (r,t), namely 

D (r, t)=E (r, t) + AttP (r, t) , where P (r, lu) ^ x (^) E (r, u) (3.1) 

with P and E the time Fourier transforms of P and E and x{uj) the fre- 
quency dependent electric susceptibility tensor. The relations H3.1|l can be 
viewed as a macroscopic picture of the interactions between the electro- 
magnetic field and the material medium, and we can naturally interpret 
the hidden variables that arise in the spectral theory of this equation as a 
model for the material medium. For the time being, however, we ignore 
such considerations (which we take up in section 5 wherein we discuss spe- 
cific models). Herein we focus on the abstract, mathematical, structure of 
the hidden variables. 

Thus given a dissipative system, defined by either the evolution equa- 
tion (|l.l|l or its frequency counterpart (|1.21|l . our main problem is to find 
a larger conservative system, governed respectively by the equation (|1.6() 
or its frequency counterpart (|1.28|l . that reduces correspondingly to Hl.lfl 
or H1.21|l upon integrating out the hidden variables. Clearly the power 
dissipation condition in the form H1.19|l for the friction function a (t) or in 
the form H2.52|l for the admittance function 21 (w) is a necessary condition 
for the existence of such an extension, as indicated by the relations H2.15t 
12.16(1 and 1(2. 5012.5111 obtained for the truncation of a conservative system. 
Remarkably, we shall see it is also a sufficient condition. 

The equality 1(2.15(1 is a possible base for the construction of such a con- 
servative extension. Thus, we pose the following problem: given a friction 
function a (t) satisfying the dissipation condition 1(1.19(1 find an operator T 
and a self-adjoint operator fii for which the equality ((2.151) holds. Then the 
desired conservative system is ((2.1() - (I2.3() . Alternatively, we may start with 
the relation ((2.5U() and ask, given the admittance 21 (C) with Re 21 {() > 0, 
whether we can find an isometric truncation T : Ti. ~^ H along with self- 
adjoint > and A for which 1(2.50(1 holds. These two closely related 
approaches each lead to constructions of an extension, based respectively 
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on operator versions of the following fundamental results: Bochner's The- 
orem, P Section 60], ^1 Section XI, 13, Theorem 2], and the Herglotz- 
Nevanhnna Theorems, P Section 59], |7|, (TUl Section 32.1, Theorem 2, 
Theorem 3]. In the second approach, the Naimark Theorem on positive 
operator valued measures plays a key role, Vol II, Appendix I, Section 
I], ^1 Appendix, Section 2, Theorem I]. 

3.1 Approach via Bochner's Theorem 

We begin by recalling, Jj Section 60], Section XI, 13, Theorem 2] and 
[Tl Theorem IX.9]: 

Theorem 3.1 (Bochner) A complex- valued continuous function s{t) of 
— oo < i < oo, is representable as 

/oo 
e-'^^dNia) (3.2) 
-OQ 

with a non-decreasing, right-continuous bounded function (a) if and only 
if s (t) is positive-definite in the following sense 



s{t-T)(p {t)tp (t) dtdr > (3.3) 



oo ^ — oo 



for every continuous function {t) with compact support. 

The measure dN (a) may be realized as the spectral measure associated 
to the vector xr and the operator fli of multiplication by a on L^(]R, dN). 
Thus 

s(t) =re-'"^ir^ (3.4) 
with r the following rank one operator 

rV- = J ^ (c^) dN (a) : L^{R, dN) C. (3.5) 

Hence, Bochner's Theorem may be seen as a special case of the following 
result, which we prove below in Section 6. 

Theorem 3.2 Let B{Hq) be the space of all bounded linear operators in 
iJo- Then a strongly continuous S (iJo)-valued function Oe (t), — oo <t< 
oo, is representable as 

ae(t) = re-'"^^r^ (3.6) 

with ill a self-adjoint operator on a Hilbert space Hi and F : iJp Hi a 
bounded linear map, if and only if {t) satisfies the dissipation condition 
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H1.19|l for every continuous Hq valued function v{t) with compact support. 
If the space Hi is minimal - in the sense that the linear span 

{gini)TK:g€Ccm, v e Hq) (3.7) 

is dense in Hi - then the triplet {_ffi,J7i,r} is determined uniquely up to 
an isomorphism. 

Remark 3.3 In fact, it is sufficient to assume that Oe [t) is locally bounded 
and strongly measurable, strong continuity then follows from (|3.6|l . 

Remark 3.4 (^g [Vli) V^v : g ^ Cc (K) , v G -ffo) denotes the linear span, 
i.e., the subspace of linear combinations of finitely many elements of the 
form g (r^i) T^v. 

Our proof of this theorem is a very elementary generalization of the 
proof of Bochner's Theorem given in ^3]. Nonetheless we are not aware 
that Theorem 13.21 has appeared previously in the literature. 

Theorem 13 . 21 provides a basis for constructing a conservative extension 
of a dispersive system without instantaneous friction: given a system de- 
scribed by a vector u in a Hilbert space Hq and governed by a dissipative 
evolution H1.6() with a strongly continuous friction function, we simply rep- 
resent it as the restriction of a conservative system in the block-matrix 
form H2.6|l with fii and F the operators obtained from Theorem 13. 21 

The construction afforded by Theorem 13. 21 is sufficiently general to de- 
scribe most systems of interest, excepting those with instantaneous friction. 
For such a system we must admit a friction "function" a (t) which is not 
strongly continuous. For most cases of interest though it is sufficient to 
assume the following condition 

Condition 3.5 (friction function) The friction function a (t) is of the 
form 

a (t) = aooS (t) + a (t) (3.8) 

where ttoo is a bounded non-negative operator in Hq and a (t) a strongly 
continuous and bounded B (_ffo)-valued function for t > 0, i.e. 

< aoo < CIho, C <oo; sup\\a{t)\\ < oo. (3.9) 

The extension (t) = 20^06 (t) + ae (t) , — 00 < i < 00 of the function a {t) 
is defined by the formula H1.18|l . 

The following result for a friction function satisfying Condition 13 . 51 will 
be useful. We give its proof below in Section 6. 
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Theorem 3.6 Suppose that the friction function a [t) satisfies Condition 
13.51 Then its extension Oe {€) = 2aoo<^ {t) + oteit) , ^oo < t < cxd, is repre- 
sent able as 

ae {t) = DUmre-""^ {T^iy , ^ (^|| + /^^^ , (3.IO) 

with Qi a seh-adjoint operator on Hi and F : I? (f^i) ^ ^0 an f2i-bounded 
hnear map, if and only if (t) satisfies the dissipation condition (|1.19|l 
for every continuous Hq valued function v(t) with compact support. If the 
space Hi is minimal - in the sense that 

((rg v. gecm, veHo) (3.11) 

is dense in Hi - then the triplet {_ffi,f2i,r} is determined uniquely up to 
an isomorphism. 

Remark 3.7 Here Dlim indicates the distributional limit, i.e. 

/C30 poo 
ae{t)v{t)dt^ lim / Te-''^^^^\T^v{t)dt (3.12) 

for every smooth Hq valued function v [t) with compact support. In essence 
we have (t) = re~'*^^r^, but this expression is ambiguous so we intro- 
duce a sort of principle value by regularizing with <I>^. 



3.2 Approach via Herglotz-Nevanlinna Theorems 

The approach via Bochner's Theorem just outlined is quite straightforward 
and adequate for many purposes, however there is an equally useful method 
which works in the frequency domain and makes use of analytic function 
theory. This second approach is based on an alternative description of the 
friction function a (t) through its Laplace transform d (C) as defined above 
in eq. 

Condition 13 . 51 for the friction function readily implies that 

a(C)=«oo+a(C), l|a(C)ll < ImC>0. (3.13) 

One advantage of d (C) over a (t) is that it is an analytic function, even if 
a (t) has the singular term aoo^ (t)- Formally the Laplace Transform d {(,) 
becomes the Fourier Transform d (w) for ^ = w with real lu. 

The power dissipation condition H1.19|l implies that Im d{() > for 
ImC, > 0, as we have seen in (|2.42f) . Therefore for each v G H, ( 1-^ 
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{v, ia (C) v) is an analytic map of the upper half plane into itself. There 
is a classical representation theory for such maps, Section 59], [7], |E1 
Section 32.1, Theorems 2, 3], which provides a tool for the construction of 
a conservative extension governed by (|2.1U|) and H2.11|l . 

Theorem 3.8 (Nevanlinna) Every analytic function g {() in the upper 
half-plane Im C > whose imaginary part is everywhere non-negative and 
which satisfies the growth condition 

limsupyy |g (ir;)! < oo (3.14) 

77 — >+oo 

can be expressed uniquely in the form 

»(o^r^, (3.15, 

where N (a) is a non-decreasing, right-continuous, bounded function such 
that 

/oo 
dN {a) = limsupTylmj^ (177)} < 00. (3.16) 
-00 ?7^+oo 

In fact. Theorem 13. 81 is a special case of the following result. 

Theorem 3.9 Every analytic function g {() in the half-plane {Im^ > 0} 
whose imaginary part is everywhere non-negative can be expressed uniquely 
in the form 

9 (C) = e + PC + r dN {a) (3.17) 

where N (a) is a non-decreasing, right-continuous, bounded function, ^ is 
real, and p > 0. 

Remark 3.10 The measures from Theorem 13.91 and Theorem 13.81 are re- 
lated by the identity dN (a) = (l + a'^^ dN (a) , and the former result is 
obtained from the later by noting that (|3.14|l implies that p = and that 
(1 + (7^) dN (cr) is a finite measure. 

Returning to the analysis of a (t), suppose first that the instantaneous 
friction vanishes, aoo = 0. Then in view of eq. H3.13|l we have 

a(C)=a(C), and ||a (C)|| < (3.18) 

Hence, given v G Hq, the function {v,id{()v) satisfies the hypotheses of 
Theorem 13.81 and consequently there is a finite Borel measure diV„_i, (cr) 
such that 

i<MOv)= r'J^. (3.19) 
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For each pair v,w £ Ho, we define the "off-diagonal" measures dN^^^ (a) 
via polarization, 

dNy^m (a) — - [dNy+tu^y+w (cr) — dNy^yj^y^yj (cr) 

(cr) + \dNy (a)], (3.20) 

so that 

{v, id (C) w) = r 'BjI:^^ for aU v,w Ho. (3.21) 

J-oo cr-C 

Using the measures dN^ ^j we define, for each ct G M, a non-negative 
quadratic form 

Qa{v,w) = f dNy^y,{a'). (3.22) 

Because 

/oo 
dNy^y (a) 
-oo 

= lim sup 77 Im {(u, ia (iyy) u)} < I sup ||q; {t)\\ j ||w||^ , (3.23) 

I'j — ^oo \ i / 

we see that these forms are bounded. Thus for each cr G M there is a non- 
negative bounded operator K {a) which satisfies Qa (v, w) — (w, K (cr) w). 
It is easy to see that this generalized spectral family of operators satisfies 
the following condition. 

Condition 3.11 (generalized spectral family) 

1. K (cr), CT G K are bounded non- negative operators in Hq- 

2. K{a) <K (A) for a < \; K [a + Q) ^ K (cr); stlim<,_,_oo K (cr) = 0; 
K (+oo) = stlimo-^oo K (cr) exists and is bounded. 

There is a fundamental result due to Naimark, ^2 in Russian], ^ Vol. 
II, Appendix I, Section I], |13| Appendix, Section 2, Theorem I, includes 
uniqueness] , which provides a canonical representation for generalized spec- 
tral families. 

Theorem 3.12 (Naimark) Let K{<t), cr g M be a generalized spectral 
family satisfying Condition 13.111 Then there exist a Hilbert space Hi, a 
bounded map V : Hi Ho, and a resolution of the identity E (cr) , cr G M 
of Hi such that: 

is: (cr) = (cr) cr G R. (3.24) 
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If the space Hi is minimal - in the sense that 

{E (a) r^u : CT e M, ve Hq) (3.25) 

is dense in Hi - then the triplet {Hi, {E (a) , a £ M.} ,T} is determined 
uniquely up to an isomorphism. 

The construction of a conservative extension from these results proceeds 
as follows. Given a (C) with aoo — 0, we obtain K (a) from Theorem 13.81 
and thence F and E (cr) from Thcorcm l3.12l Letting fli be the self-adjoint 
operator adE (cr), we obtain 

poo -, 

a{Q) = -\ -TdE{a)T^ 

J-oo cr - C 

/OO 
adE (a), (3.26) 
-oo 

which is the desired representation for a (C) ~ compare with eq. H2.54|l . 
This construction is summarized by the following operator generalization 
of the Nevanlinna Theorem. 

Theorem 3.13 Every S (ifo)-valued analytic function G {() of the upper 
half plane Im^ > with ImG (<^) everywhere a non-negative operator, and 
which obeys the growth condition 

limsup77|jG(ir;)|| < oo, (3.27) 

7] — *4-oo 

can be expressed in the form 

G(C) = r[f7i-affJ^'rt (3.28) 

with fli a self adjoint operator on a Hilbert space Hi and T : Hi ^ Ho a 
bounded map such that 

TT^v = lim rfG (ir/) v for every v € Hq. (3.29) 

If the space Hi is minimal - in the sense that 

(/ (Qi) T^v. / e Cc (R), ve Ho) (3.30) 

is dense in Hi - then the triplet {Hi,D,i,T} is determined uniquely up to 
an isomorphism. 

Remark 3.14 Theorem l3 . 1 3l and operator generalizations of the full Herglotz- 
Nevanlinna Theorem 13 . 91 are certainly known to experts, and have seen ap- 
plication, for example, in the theory of self-adjoint extensions of symmetric 
operators - see, e.g., @]. 
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Remark 3.15 In the scalar case, the function (/(C) in (|3.15|) may be ex- 
pressed in terms of the resolvent of the self-adjoint operator fii^ (a) — 
(J(I>{cf) on {dN), 

C/ffJ^'rt, (3.31) 

T : L^{R,dN) C (3.32) 

the rank one operator as in the discussion of Bochner's Theorem above. 
Thus Theorem 13.81 may be seen as a special case of Theorem 13.131 (Al- 
though the latter is in fact a consequence of the former.) 

For systems with instantaneous friction (aoo 7^ 0) , we shall make use of 
the following result, which is derived from Theorem 13 . 91 and the Naimark 
representation fTheorcm l3.12|l in the same way as Theorem 13. 131 We dis- 
cuss the proof of this result, and related constructions, in Section 7 below. 



with 



g(C)-r[r!i- 

= / V (cr) dN (cr) , 



Theorem 3.16 Let G (C) he a B (iJo)-valued analytic function of the up- 
per half plane ImC > with ImG (C) > everywhere. If 

lim G (ii?) V ~ igooV (3.33) 

R — >oo 

exists for all v G Hq with goo > a non-negative operator, then G(C) can 
be expressed in the form 



GiOv - lim r [n, - CIhA ' [ ) v, for aU V e H, (3.34) 

with ill a self-adjoint operator on Hi, and F : V (ili) — > Hq a fii-bounded 
linear map. If the space Hi is minimal - in the sense that 

((F/ (ni))^ v. f eGciR), ve Ho) (3.35) 

is dense in Hi - then the triplet {i/i,i7i,F} is determined uniquely up to 
an isomorphism. 

For a friction function a (t) satisfying Condition 13 . 51 and . consequently, 
a (C) satisfying (|3.13|) the fimction G {Q = id {() satisfies the hypotheses of 
Theorem l3.16l with 500 = Q^oo- Hence, we can apply Theorem l3.16l to obtain 
operators F and f2i which provide a basis for a conservative extension via 
the abstract model of the previous section. The restriction that ^oo be self- 
adjoint is natural in this context, since any anti-Hermitian contribution 
to lim/{_>oo d (ii?) can be absorbed in —iA. In the present context, the 
operator G appearing in Prop. 12.21 is G — d{i), which is in fact bounded. 
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3.3 Spectral representation of the admittance opera- 
tor 

The admittance operator 21(C) = i [C"^ — A + ia associated to the 

dispersive system (jl.lf) is an analytic function on the domain {Im^ > 0}, 
and the power dissipation condition takes the form Re 21 (C) > when 
expressed in terms of 21 (see eq. (12.5211 .') Furthermore, under condition 
13.111 we see that 

hm TySt (ir;) = TO^^. (3.36) 

?7 — *oo 

Thus i2l (C) satisfies the conditions of Theorem 13 . 1 31 and therefore has the 
following representation 

2t(C) = ir(c-f7)"^r^ (3.37) 

with n self-adjoint on a Hilbert space H and T : H ^ Hq bounded. In 
this section we develop a basis for a slightly different representation of the 
admittance operator, namely 

21(C) =iT(CX-yi)"'r^ 

with T : Ti ^ Hq an isometric truncation. (3.38) 

In fact, eq. (|3.38|) is a consequence of eq. (|3.37|) and the polar decom- 
position r = TAl"5where T : H ^ Hq is an isometric truncation and 
< 6 < M. The relationship between eq. (|3.38|l and eq. H3.37|l is 
expressed through the identity A — Ai^Q,/^! . However, there is some 
flexibility in the choice of M, which we discuss below. Since 
may not be well defined for unbounded M, we shall require throughout 
that M be bounded. 

Let us first fix some notation. Given two Hilbert spaces H, Hq and a 
bounded linear operator L : 7i — > Hq we denote 

KeiL ^ {V en : LV ^0} , Rani = {LV -.V eH}. (3.39) 

We denote the closure of a subset 5* in a Hilbert space by S, and the 
restriction of L to 5* by L\g. We need the following elementary facts 

KerL = [Rani^]^ , RalTX = [KcrL^]^ , (3.40) 

where [•]"'" denotes the orthogonal complement in the relevant Hilbert space. 
We refer to the orthogonal direct sum decomposition 

n = RanLt e KerL where Ran it ^ [KerL]"^ , (3.41) 

as the L-decomposition of the Hilbert space TC. 
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Theoreni l3 .131 readily implies the following statements regarding the de- 
composition associated to an operator function G (C) of the type considered 
there. 

Corollary 3.17 Let G{() : Hq Hq satisfy all the conditions of Theorem 
IXni so G (C) = r (r^i - (Iny^ rt. Then we have 

KerG(C) = Ker G^(C) = KerT^ for all ImC > 0. (3.42) 

Furthermore, the F^^-decomposition of the Hilbert space Hq, i.e. 

Ho = Ho® Ho , Ho= [Ker F^] ^ Ra^H^ (3.43) 
reduces G(C) in the sense that G(C) has the following block form under 

G(C)-G^„(C)©0= [ I , 

where G^^(C) = f [fii - (/ffj"' : Ho ^ Ho, f = P^T. (3.44) 
In addition 

KerG^^(C) = Kerf^ = {0} , Ranf = Ho, ImC > 0. (3.45) 

Corollary 13 . 1 71 allows us to extract from the operator function G{() its 
nontriyial component and motiyates the following definition. 

Definition 3.18 For any function G(C) : Ho Ho satisfying all the con- 
ditions of Theorem 13 . 1 31 the pair ^Ho,Gfj^^{C)^ defined in Corollary 13.171 
is called its reduced representation. 

For the representation eq. (|3.38(l to hold we shall of course require that the 
function G (C) = — iSl (C) is in reduced form. Howeyer, to obtain a bounded 
mass operator, we must strengthen this requirement by assuming there is 
e > such that 

ImG(i7?) >-,?]> I, (3.46) 
V 

which implies that FF^^ > e. Under this additional assumption, we haye 
the following theorem, which proyides a basis for a spectral representation 
of the admittance operator. 

Theorem 3.19 Let G (C) be a S (i7o)-yalucd analytic function of the up- 
per half plane Im C > with Im G (C) > for eyery C and assume that 
G (C) obeys the growth conditions eq. (|3.4t)|) and 

limsup77||G(i77)|| < oo. (3.47) 
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Then there exists a Hilbert space H, such that 

G(C) = r [17 - (^InV^ where is self-adjoint in 

T Hois bounded, and TT^ = Hm rjG {it]) > e. (3.48) 

rj — ^+00 

Furthermore, if T denotes the isometric truncation T ;= [FF^] ^ F, then 
T = UPh'^ with U = T|^/ a unitary map from H'^ — [KerF]^ to Hq, and 
the following representation holds: 

G{C)=T[Am~CM]-^T^ , (3.49) 

where Am = M^VlM^ is self adjoint in H and M : H H is any 
operator of the form 

M ^mG®mi, (3.50) 

with mi : KerF ^ KerF bounded and strictly positive (mi > SIkbtT, 
6 > 0) and 

ma = U-^ [rr^y^ U. (3.51) 

Note that is bounded and strictly positive since mg > llFFtjl 2 
and llmdl < j-- If the space TC is minimal - in the sense that 

{fmr^v: /eC,(R), veHo) (3.52) 

is dense in 7i ~ then the triplet {7Y,A^,^} is determined uniquely up to 
an isomorphism and the choice of mi. 

As the reader may easily verify. Theorem 13.191 follows from Theorem 
13.131 and the polar decomposition (see, e.g., ref. j6j Section VI. 7]), summa- 
rized here: 

Theorem 3.20 Let T : H ^ Hq be a bounded linear operator with 
KerFt = {0} and let 

n = H'^QHi, [KerF]^ , Hi = KerF (3.53) 

be the F-decomposition of 7i. Then 

TV for any nonzero V H'^ and VW^ = Hq, (3.54) 

F^F = KPh' where K = F^F| > 0, (3.55) 
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and the following "polar decomposition" holds 

r = UK^Ph'^ where U -.H'q^Hqis unitary, (3.56) 

indicating, in particular, that H'q is an isometric copy of Hq. In addition, 

rrt = uKU-\ K = r^TL, = u-^tt'^u, (3.57) 

T=[VT^YUPh',, (3.58) 
and r is an isometric truncation if and only ii K = Ih'^ , in which case 

r = UPh'^ . (3.59) 
Furthermore, let K be any operator in H of the form 

lC = K®Ki with Ki:Hi^ Hi non-negative. (3.60) 
Then the following "generalized polar decomposition" holds 



r = T/C^ = T 



K2 ® 



(3.61) 



where T = UPr'^ ■ ^ Ho is an isometric truncation. Notice that T and 
K are uniquely determined by F, but /C depends on the choice of Ki. 

4 Final schemes of the construction of the 
conservative extension. 

In this section we summarize in a concise form the main consequences 
for time dispersive systems of the analysis of preceding sections. There 
are two intimately related ways to construct a conservative extension of a 
dispersive system. The first, the "friction function scheme," is based on 
the operator- valued friction function a{t) or its Laplace transform a{C,), 
whereas the second, the "admittance operator scheme," is based on the 
admittance operator 21 {(). 

4.1 Friction function scheme 

The original time dispersive system is described by an evolution equation 

mdtv (i) = -iAv (i) - / a (r) v {t - t) dr + f (t) , v {t) e Hq, (4.1) 
Jo 

with Ho the space of system states v describing "observable" variables, A a 
self-adjoint operator in Ho describing the internal dynamics, a (r) : Hq — > 
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Hq an operator-valued friction function accounting for time dispersion and 
losses, and / (t) £ Hq a time dependent external force. Based on the 
notion that the friction function a (r) is a result of a coupling between the 
"observable" variables v Q Hq and some "hidden" variables described by a 
vector w belonging to a Hilbert space Hi , we seek a conservative ( not time 
dispersive) extension of eq. i4-l\ ) in the form 



mdtv (t) = -iAv (t) - iTw (t) + f (t) , 
dfw (t) = -iT^v (t) - ifiiw (t) , 



(4.2) 
(4.3) 



where fli is a self-adjoint operator in Hi describing the internal dynamics of 
the "hidden" variables and F : i/i — > Hq is a coupling operator between the 
"hidden" and "observable" variables. The system (|4.2|I - H4.3|I . of course, is 
expected to reduce to the original system (|4.1|) after the "hidden" variables 
w (t) are eliminated by solving eq. (|4.3I) and plugging the result into eq. 
(|32J|. 

Now, the question is how to construct the conservative extension of the 
form H4.2|l - H4.3|l based on the original dispersive system (|4.1|l ? As we have 
seen, a necessary and sufficient condition for such an extension to exist is 
the power dissipation condition 



{v (t) , {t — t) V (r)) dtdr > 0, for every v{t). 



(4.4) 



— OO ^ — OO 



A first possibility is to obtain the triplet {Hi,r, fti} from the function a (t) 
via the operator version of the Bochner Theorem l3.2[ which represents a (t) 
as 

a(t) =re-""ir^ (4.5) 

However, in practice most systems are specified in the frequency domain, 
and it is more convenient to carry out the construction in that setting. 

Taking the Fourier-Laplace transforms of (|4.1|) and (|4.2|) - (|4.3|) with re- 
spect to t yields 

mCv (C) = [A~ia (C)] V (C) -f i/ (C) , C = ^ + i^, = Im C > 0, (4.6) 
and the conservative system, recast in block matrix form, 

(^MV (C) = AV (C) +iF{C)., C^uj + iry, r; = Im C > 0, (4.7) 



A = 



A r 
rt fii 



, M = 



wiC) 

m 
Ih^ 



/(C) 





, v{0=TV{0, T=[ Ih, ] 
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Note that the operator T defined in H4.7II is an isometric truncation of the 
extended space ® Hi onto the space of observable variables i?o- After 
the Laplace transform, the power dissipation condition (|4.4|) becomes 

Rea(C) > for ImC > 0, (4.8) 

which is equivalent to (|4.4() . 

The scheme of the construction is as follows. Given an operator-valued 
friction function d (C,) which satisfies the power dissipation condition H4.8|l . 
we apply the operator version of the Nevanlinna Theorem formulated in 
Theorems 13.131 and 13 . 1 61 to construct a triplet {iJi,r,f2i} giving the rep- 
resentation 

a(C) = ir(c/H, -i7i)-'rt, (4.9) 

T : Hi ^ Hq, fli : Hi ^ Hi is self-adjoint, 



ii\\a{irj)\\^0{rj-^), 



a(c) = stiimir(c/H.-f^O-^^r^^^-^j , (4.io) 

if stlim,,^oo (i*?) exists. 

In more detail, assuming \\a (i»y)|| = 0{ri^^), we first obtain the operator 
version of the Nevanlinna representation 



a(C)-i/ 7 dK{a), (4.11) 

where dK (cr) : Hq Hq is an non-negative operator-valued measure over 
R. Notice that dK (a) is not necessarily a resolution of identity, and for 
different intervals Ai and A2 in M the operators K (Ai) > and K (A2) > 
may not commute. However, having found dK (a) we apply the Naimark 
Theorem 13.121 and get (i) a Hilbert space Hi ; (ii) a resolution of identity 
dE (ct) ; (iii) an operator T : Hi ^ Ho such that 

dK{a) ^rdE{a)r^. (4.12) 

To complete the construction of the triplet {iJi,r, fii}, we define 

/C30 
adE {a). (4.13) 
-C30 

It is when applying the Naimark Theorem \3.1'M hat we get the desired triplet 
{i?i,r, ill}, which is the central point of the construction of a conserva- 
tive extension for a time dispersive system satisfying the power dissipation 
condition \4-4\) , \4-^ - 
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We may present a more constructive picture of the triplet the triplet 
{Hi,T,fli} under the additional assumptions that 

Re a (o- + iO) = stlim Re a (cr + i??) (4-14) 

exists for almost every a, and 

lim / ||Rea(CT + i7/)-Rea(cr + iO)||dcr = (4.15) 
v^o J J 

for any finite interval /. With these assumptions, the Stieltjes-Inversion 
formula - discussed in Appendix A.l below - provides the following explicit 
formula for dK: 

dK {(j) = i Re a (cr + iO) da. (4.16) 

TT 

Let us define N (a) — tt^^ Re a (cr + iO), and note that N (a) is a non- 
negative operator for almost every a. We choose the Hilbert space Hi = 

(K, Hq), which is the space of square integrable functions from M to 
with inner product 



{wi,W2)h, = / {wi (a) , W2 {(^))ho (4-17) 

J — oo 

Choosing for Q.i the operator 

n^w [a] = aw {a) , (4.18) 

and for T the map 



oo 



Tw= j y/N {a)w [a) da, (4.19) 

J — oo 

it is easy to see that 

pea -, 1 

a(C)=i/ N{a)da = \T——T\ (4.20) 

which is the desired representation. 

The above representation on Hi ~ I? (K, Hq) may not be minimal if, 
for a from a set of positive measure, A'^ (a) has a non-trivial kernel. In that 
case it is more useful to consider the Hilbert space Hi ~ I? (dK) defined 
to be the space 

/>oo ^ 

{w (a), N (a) w (a)) da <oo} , (4.21) 



w : M ^ Hq measurable 

modulo null functions with A^ (a) w (a) ~ for almost every a. In fact the 
proof of the Naimark theorem (without assumptions l|4.14|l - 14.15|l ) proceeds 
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by constructing the space LP' [dK). A sketch of this construction is given 
in Appendix A. 2. 

Finally, having found the triplet {iJi,r, f^i}, we construct the conser- 
vative extension as the system of equations (I4.2|) - (|4.3|) . Using the matrix 
form 14. 7|) of the conservative system (|4.2|I - H4.3|I we also get the following 
representation for the admittance operator 

21 (C) = iT [CM - Ay^ = {Cm -[A-ia (C)]}"' . (4.22) 

The above construction oj the triplet {_ffi,r,f2i}, including the Hilbert 
space Hi of "hidden variables ", is essentially independent of the operator 
A which describes the internal dynamics of the "observable" variables. 

4.2 Admittance operator scheme 

Suppose the original time dispersive system is given by its admittance 
operator 2t(C), acting in the Hilbert space of "observable" variables Hq 
by the following equation 

(C) = 21 (C) / (0 , C = ^ + ri = Im C > 0, (4.23) 

relating the generalized velocity v {() and the generalized force / {() in the 
complex frequency domain. For real ( — u; the equation (|4.23() reduces to 
the familiar real frequency form 

v{uj)=^{Lo)f{Lo). (4.24) 

We assume that the Hilbert space Hq and the admittance operator 21 (C) 
are already reduced in the sense of Definition 13.181 i.e. Hq = Hq and 
2tj^ (C) = 21(C), that 21(C) satisfies the power dissipation and the 
growth conditions 

Re2t(C) > for ImC > 0, limsupyy |j2l (i?7)|| < oo. (4.25) 

77 — >+oo 

We seek a conservative, non time- dispersive extension, of eq. \4-^'i^ via 
the following representation for the admittance operator 

^{0 =iT{CM- Ay^T\ (4.26) 

with T an isometric truncation from a Hilbert space Ti to Hq, A4 > SI-h, 
6 > on H, and A self-adjoint in H, corresponding to the following evolu- 
tion equation for v, 

MV {t) = -lAV (i) + tV (i) , v{t)= TV {t) . (4.27) 
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To construct the representation (|4.2()|) , we first construct the Hilbert space 
7i and the representation 

21(C) = i^a(C2:-^^a)"'4, (4.28) 

following the argument which led to 14. 9|) for the friction function in the 
previous section - applying Theorem l3.19l to G{() — i2t(C). 

Having obtained (|4.28l) . we use Theorem l3.20l to get the following polar 
decomposition 



Ta = Um^'^PH'^ where U : ^ Hq in unitary, (4.29) 



ma = U 



-1 



-1 



u > 



-1 



Im, ra4= lim iTySl (iry) > 0. 



Then we introduce the mass operator 

M = ma © mi where nii : Hi Hi > SIhi , S > 0, (4.30) 
with nil being chosen as we please. In view of (I4.29f) 

Ta = UPhiM--- = UPhi^ (m^^' m^') , (4.31) 
and we get the desired representation 

21(C) = vr{QM-Ay^T\ 

where T ^ UPh'^ said A = U^^M^ . (4.32) 

Note that the Hilbert space of "hidden variables" is Hi = H Q H'q, which 
is kerFa, i.e. 

Hi^iijeH: Ta^A - 0} . (4.33) 

5 Examples of the construction of conserva- 
tive systems 

In this section we apply the general scheme for the construction of conser- 
vative extensions, as described in Sections 3 and 4, to a few well known 
classical systems: a damped oscillator, a general scalar dispersive dissi- 
pative system, and a classical dielectric medium. Related constructions 
have appeared elsewhere in the literature, e.g.. Lamb's representation of a 
damped oscillator as a mass attached to an infinitely long tense string [5]. 
Interesting examples and very detailed studies of relations between admit- 
tance operators and spectral measures for loaded strings (as described by 
Krein- Feller operators) are offered in the second paper in ref. [Jj, wherein 
dispersion is introduced via boundary conditions. 
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5.1 Damped oscillator 



We consider a damped oscillator with mass ttIq > 0, real frequency fio, and 
friction coefficient 7o > 0, described by a complex variable v which evolves 
according to the following equation 

niodtv ^ -imofiov - ^oV + f (t) , (5.1) 

where f (t), the external force, is a complex- valued function. Evidently, 
the equation 1)5. l|l is a particular case of a system of the form with 



Ho = C, A = nio^o, 

a (t) = Qo (t) = 7o<5 (t) , and do (C) = 7o for ImC > 0. (5.2) 

Following the friction function scheme described in Section 3, we look 
at eq. (|5.1() as a consequence of a larger conservative system with hid- 
den degrees of freedom of the form (|4.2|l - (|4.3(l . and proceed with the con- 
struction of the triplet {Hi,ili,T} based on Theorem 13.161 Notice that 
G {() = iSo (C) — 170 evidently satisfies the hypothesis of Theorem 13.161 
Hence, we must have 



7o = i hm r— — 



i?2 



t 



, ImC > 0. 



(5.3) 



The discussion following eq. H4.15|l - in particular eq. (|4.2(J|I with 

N (cr) — 'Jo/tt - suggests that we take Hi = (R), 



rJiV- (cr) =a^p{a), cr e M, V (cr) e L2 (K) 



(5.4) 



and 



— lim 



:i>(<^)> 



where 1 = 1 (a) 



1, cr e 



(5.5) 



Although r is not bounded, we note that F is an fJi-boundcd map from 
V (ill) ^ C, since 



|F^|<,/^^lim 





1 


J-R 


a + i 



\ia + i),Pia)\da<^ || (fli + i) V'll • (5.6) 



Formally, 



(rt„) (a) = y^t;, a e 



(5.7) 
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The generator of the dynamics in the conservative system is the self- 
adjoint operator A acting in the Hilbert space H = CS) L2 (K, da), defined 

by 



A 



V 



2^1 



V 



mo^oV + (1, ^ (cr)) 

^/^v + aijj (cr) 



(5.8) 



cr e 



with the domain 
D{A) = 



V 



cr e 



V + a^ (cr) e L2 (M, cfcr) 



(5.9) 



One can show that for any vector in D [A) the expression (1, ^ (cr)) is well 
defined as the following limit 



lim 

00 



(cr) dcr, 



(5.10) 



and by Prop. 12.21 we see that A is self adjoint. 
Finally the mass operator here takes the form 



M = 



nio 
Ih, 



(5.11) 



and the desired extended conservative system for the damped oscillator is 



ruodtv 



brio^oV — U — / ip (a) da + f (t) 



5tV i'T) = 
for iv,Tp) e n = C(SL2 (M, da) 



iatp (cr) + i\l 



(5.12) 



5.2 General scalar dispersive dissipative system 

A number of classical dispersive systems, including homogeneous dielectrics 
and an oscillator with a retarded friction, can be described by a scalar 
complex variable v governed by an evolution equation 



rrisdtv (t) = —imsQsV (t) 



a, (t) v{t-T)dT + f (t) , 



(5.13) 



with fls (r) complex-valued functions, rris a positive number and fls a real 
number. The function / is the external force. The scalar friction func- 
tion flg (t) is assumed to satisfy Condition 13.51 and the power dissipation 



36 



condition (|1.19|l . which is to say it is a positive definite function as in the 
classical Bochner's theorem. The system described (|5.13|) has, according 
to (|2.44|) . the following admittance form 

^'(0 = a, (0/(0, 

2ts(0 ={TOs(C-aO + ia4C)}"\ ImC>0. (5.14) 

Let us consider here a general scalar dispersive dissipative system given 
in admittance operator form 

V (0 = 2ls (0 / (0 , C = ^ + iv, V = Im C > 0, (5.15) 

where the scalar admittance operator Slg (0 satisfies the power dissipation 
condition 

Reas (0 > 0, C = w + i77, 77 = ImC > 0, (5.16) 

and 

lim sup r] |2ts (i?7)| < 00, (5-17) 

and 2ls (C) 7^ at least for one 

To find a conservative extension for (|5.15|l . we use the admittance oper- 
ator scheme from Section 4. In this case Hq = C and we first seek a Hilbert 
space H and a self-adjoint operator fia^ in it such that 



2ls (0=ira. (a-f^mj"'4. 



r°° Ta En,, (da) tJ, 
/ ° V ; g, ^ im^>o. (5.18) 

J -00 CT-C 



We construct the spectral representation H5.18(l as follows. In view of con- 
ditions (|5.16(l and H5.17(l , the classical Nevanlinna Theorem 13.81 gives a 
non-negative scalar measure dN<2i^ (a) satisfying 

i2t, (0 = r > 0, dTVa, (a) > 0. (5.19) 

Let us define by 
1 



dA^a, {<y) = Jim t^^s (ir]) , (5.20) 

and note that if 2ls is given by H5.14(l . ma^ — mg. 

We set 

n^L2 (dA^a^ (o-) , C) , H^^{Tpen:tpia)=veC, aeR}. (5.21) 
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In other words, H'q is defined as a set of constant functions of a, which 
is, evidently, unitarily equivalent to the set of complex numbers C through 
the following mapping 

U : H'n-^C, U (vl) ^ , where 1 = 1 (a) = 1 for e R. (5.22) 

We also define 



[Fa,^] = (1, ^) = (a) diVa, (a) e C, (5.23) 

implying 

r^^w = wl, ve C. (5.24) 

The representation (|5.18() readily follows from the definitions H5.21|l - 
Also, from if?:^ and if^?^ it follows that 



a. 



dN<^^ (a) 



Observe that the orthogonal projection Ph' is 



m<2L^ / V (ct) (iA^a, (c) 



1, ere 



(5.25) 



(5.26) 



We may choose the mass operator M to be just the scalar operator M = 
rriQi^lT-c, and, consequently, get the desired components of an extended con- 
servative system 

H = L2 (diVsi, (fi) , C) , m^Jn, 

{a) = m^ai^ (a) , T ^ U Phi, , (5.27) 

where U is defined by H5.22|l . 

The measure dN<^ (cr) can be recovered from Stieltjes' formula - see eq. 
I|A.3|I below - 

/ (a) dA^a. (ct) 

= lim - / / (a) Re {21, (a + i??)} da, for / e Co (M) . (5.28) 

In particular, if dN%^ (a) has a density n<^^ (cr) with respect to the Lebesgue 
measure, i.e. 

dA^a^ (cr) = na, (cr) da, (5.29) 
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then the density (a) is the following pointwise limit 



"•a, (c) = lim - Re {2ts (a + i??)} , 

1)^0 TT 



(5.30) 



which, in view of (|5.2U|) . satisfies 

.00 

(it) > and / na (cr) dcr = lim yySts (i?]) = . (5.31) 

J-oc 'Tia, 



When the measure dN<^^ (a) has a density ria, (cr), there is a slightly 
different description based on the Hilbert space TC = L2 (M), namely: 



H, 



n<^^t^ (a) ^ ai, (a) , ^ e L2 (R) ; 
= {V^ e L2 (C) : ^ (d) = i;v/"2i. (^), ^' e c} 



(5.32) 
(5.33) 



Pa.V'] (^) 
Then we get 



(a) = w^na, (cj). (5.34) 



[^^//iV'] (-) = 



V"-a, (o-), cr e 



and, consequently, 

H = L2 (C) , X = majw, W) = ma,a^ (ff) , ^ = C/P^ 



(5.35) 
(5.36) 

(5.37) 



Finally, the evolution of the extended conservative system, as described 
by its state '0 (i, cr) , is governed by the following equations 

m^^dtil^ {t, a) 

= -ima^crV' {t, cr) + / (0 v^"^a,»^a, (o-), -0 (i, cr) e £2 (C) , (5.38) 

with na, (cr) = hm^^o ^ Re {2ts (a + i??)} > and / (t) the external force. 
The state v (t) of the scalar dispersive system f5.13\l is represented by the 
following linear functional of (t, cr) 



vit)^T4'{t) 



V"^a,"as (cr)V' (t, cr) dcr. 



(5.39) 



It is of interest to note that under this construction, the extended sys- 
tem described by -0 is always governed by the canonical evolution equation 
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H5.38|l in the Hilbert space H. ^ L2 (C) with generator Atp (a) = aip (c) . 
Consequently, the only feature which distinguishes different scalar disper- 
sive systems is the mass rn%^ and the "observable variable" v (t) obtained by 
projecting onto a one-dimensional Hilbert space Hq spanned by the vector 
A/mgi^ria^ (tj) in L2 (C) . Observe that the external force / (t) \/m%^n<}i^ (a) 
is in the space H'^. 

5.3 Maxwell equations for lossy and dispersive media 

In this section we construct a conservative extension of the Maxwell equa- 
tions for a homogeneous, lossy and dispersive medium following the friction- 
admittance scheme from Section 3.3. 

The classical Maxwell equations for a homogeneous, lossy and dispersive 
medium are Pl Section 1.1] 

V X E (r, t) = -~dtB (r, t) - AhIb (r, t) , V • H (r, i) = , (5.40) 



V X H (r, i) dtT) (r, t) + AnJ (r, t) , V • D (r, t) = , (5.41) 

where H, E, B and D are respectively the magnetic and electric fields, 
and magnetic and electric inductions, and J and Jb are respectively the 
external electric and magnetic currents. For simplicity, we consider here a 
non-magnetic medium, which amounts to taking 

H = B and Jb = 0, (5.42) 

in units such that /io = 1- We also assume there are no free charges, which 
is the assumption that the current J (r, t) is divergence free 

V-J(r,i)=0. (5.43) 

The dispersive properties of the medium come through the material 
(constitutive) relations which, in the simplest case of a homogeneous and 
isotropic medium, take the form 

D(r,i) =E(r,i)+47rP (r,t), 

POO 

where P (r, t) = / x (t) E (r, i - r) dr, (5.44) 
"'0 

and x{t) is the scalar- valued response (aftereffect) function, "3*, Chapter 
2] , [SI Section 3] . In the frequency domain the relation 15.44|) between the 
polarization P (r, t) and the electric field E (r, t) becomes 

P (r, u;)^x H E (r, cu) , xH= X (0 e'"* dt, (5.45) 

Jo 
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where x (w) is the so-called frequency dependent electric susceptibility, which 
is a scalar-valued function for the case we consider. Since the medium is 
homogeneous and isotropic, x does not depend on r and we have 



V • D (r, i) V • E (r, t) = V • P (r, <) = 0. 



(5.46) 



To construct a conservative extension we begin by recasting the Maxwell 
equations for the dispersive medium in the general form from Section 3. 
The "observable" variables v and the corresponding Hilbert space Hq in 
this case are 



vit)^ 



E(r,t) 
B(r,i) 

where Hq ^ {v e 



V-E(r,i) = V-B(r,i) =0}. (5.47) 



In other words, Hq consists of square-integrable 6-dimensional fields with 
the components E (r, t) and B (r, t) divergence free. Comparing equations 
H5.4U|I - (|5.46I) with the general evolution equation Hl.lll we set: 



mm 



I Ho 



A 



M 





-iV> 



iV^ 




OM (i) = fls (t) 
where Og (t) = 'i-rrdtx (*) , h 



h 



1 

1 

1 



(5.48) 



or, in the complex frequency domain 



OM (C) = as (C) 



h 




as(C) = -4^iCx(C), IniC>0. (5.49) 



The reduced form of om {(,) - as in DcfInition l3.18l - is 

H^ = [ve (C^) ; V • E (r, i) = 0, B (r, t) ^ 0} , 
V^JC)E(r)=a, (C)E(r). 



(5.50) 



Notice that the friction function am (t) is related to the time derivative 
of the polarization dtP, represented by the scalar function (C) in the 
complex frequency domain. That is consistent with the physical fact that 
the work per unit time done by the electric field E to produce the polar- 
ization is given by dtP ■ E. In view of the simple structure of the operator 
friction function H5.49(l the power dissipation condition turns here into the 
following condition for the scalar function {() = — 47ri('x {(,) 



Reas(C) =47rIm{Cx(C)} > 0, ImC > 0. 



(5.51) 
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It is easy to see that the construction of a conservative extension is es- 
sentially reduced to the construction of the conservative extension for the 
scalar friction function as (t). Thus, a conservative extension for the elec- 
tric polarization can be found using the results of the previous section for 
a general scalar dispersive dissipative system and we obtain, in particular, 
the evolution equations H5.38|l and the representation (|5.39|) . 

For simplicity, suppose that the following limit exists for every real a 



(cr) — lim — Re {as (cr + ivy)} 

ri-*+0 TT 



lim 4Im{(cr 4- ir/) X (cr + i?7)} , (5.52) 



and that 



lim ?7 Re {os (iyy)} = lim AijlmKa + irj) x {a + irj)} < oo, (5.53) 

rj — *oo ?7 — >oo 

conditions which are satisfied in many non-trivial examples. Then the 
desired conservative extension of the original Maxwell equations ^5.40^ - 
15.^6] ) for a dispersive and dissipative dielectric medium takes the form 

aiH(r,t) = -VxE(r,i), 
5tE(r,t)=VxH(r,i) 

-/ \lmxnx (c^)* (r, t, cr) da - 47rJ (r, t) , (^-54) 



m^dt^ (r, t, a) = — im^cr* (r, t, a) + \/ m^n^ (o')E (r, t) , 
where 



n-j^{a) = —A lim Im {(cr + 177) x (cr + i??)} > 0, = / n^{a) da, 
47rCx (C) = r Im C > 0, * (r, t, a) G L2 (M; C^) , (5.55) 

and the fields H (r, i), D (r, t) are divergence free, i.e. 
V-H(r,t) = 0, 

V-E{r,t)+4TTf I ^m.sjix (cr)V • * (r, r, a) dadr = 0. (5.56) 
The electric polarization P (r, t), its time derivative dtP (r, t) and the elec- 
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trie induction D (r, t) are now defined by 

J — OO J — OQ 

poo 



dtP (r, t) = ^J ^m^r^)^ (r, t, a) da, (5.57) 

D (r, t) ^F,{r,t) + f f \J m^n^ (c)* (r, t, ct) dcrdr. 

We note that according to H5.54|l the vector ^ * (r, r, a) dr in 
H5.57|l evidently represents a time dependent microscopic dipole of mass 
m^, localized at r, which oscillates with natural frequency a. The total 
polarization P (r, t) of the medium at point r is a superposition of a (contin- 
uum) number of microscopic dipoles localized at r. The natural definition 
for the energy of the "hidden" medium is 



Iff \^{r,t,a)\' dadr, 

^ JE3 J-OQ 



(5.58) 



consistent with the evolution equation and the assumption that J^s E(r, t) ■ 
dtP{r,t)dr is the instantaneous rate of work done by the electromagnetic 
field. 

In a forthcoming paper, we shall discuss in greater detail conservative 
models for dispersion in dielectric media - including scattering theory and 
inhomogeneous media - however, to complete the discussion here it may 
be useful to consider one physically relevant non-trivial example. For this 
purpose, let us take the so-called Lorentz medium in which 17, Section 
3.5], Section 9.1] 



x(t)=XL(t)=^^xp{-|i}^, i>0; (5.59) 

x(C)-xl(C)- 2 ■ r (^-60) 

with LUp, uo and 7 positive parameters. Observe that 

Im iCXL (C)} = — ; ; ^3 TTZ ^2 ' C = ^ + (5-61) 



and hence 

Im{CXL(C)}>0, ImC = 77>0, (5.62) 
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in full compliance with the power dissipation condition H5.51II . In addition, 
from H5.30|l . (|5.31l) we see that 

2 2 

n^, (a) = 4 lim Im {(a + ir;) xl [a + iv)} = 4 ^ ^ ""^V , ^ (5.63) 

''^+0 (tjg - (t2) + 0-272 

exists pointwise, and 



"^yr ^ / "-x (^) = 1™ ^7 XL (i?7) = 47rw (5.64) 

Plugging in the above values ("') ^^^d m^^ i'^to extended Maxwell equa- 
tions (|5.54|I - H5.57|I we get the desired conservative description of the Lorentz 
medium. 



6 Dissipation and continuity of the spectrum 

It may seem startling that dissipation, i.e., losses, can arise when we trun- 
cate a unitary evolution. However, such results are familiar from the theory 
of unitary dilations of contractive semi-groups |12j. which is the theory of 
solutions to for friction without retardation, i.e., a{t) — aoaS (t). Fur- 
thermore, the mechanism at work is physically very natural: there is energy 
transport to a large number of "invisible" degrees of freedom, i.e. to "heat." 

Mathematically, a rigorous analysis of losses in very large but finite sys- 
tems is generally complicated by the fact that systems with a finite number 
of degrees of freedom eventually (perhaps after an extremely long time) re- 
turn arbitrarily close to their starting configuration (Poincare recurrence) . 
However, if we consider an idealization in which there are infinitely many 
hidden degrees of freedom, resulting in infinite recurrence time, then the 
Poincare recurrence may not occur and there is hope of describing losses 
in a cleaner and simpler way. As it turns out, a sufficient condition for 
losses is strict positivity of Kea{C), which implies absolute continuity of 
the spectral measure for the generator A of the dynamics of a conservative 
extension. 

To state a quite general condition, we use the notion of non-tangential 
boundedness. Given wq S M let the cone of aperture 6 e (0,7r) at wq, 
denoted Tg{u}o), be the set 

Tg{ujo) = |w + i?7 : ?7 < 1 and \uj — ujo\ < 77 tan - j- . (6.1) 

A function F : C+ C, with C+ = {C : ImC > 0}, is said to be non- 
tangentially bounded at luq (zM. if 

sup |F(w)|<oo, (6.2) 



44 



for some 9 G (0,7r). We shall use the following well known theorem from 
harmonic analysis regarding boundary values of non-tangentially bounded 
harmonic functions, see e.g. 

Theorem 6.1 Let F : C+ ^ C be a harmonic function and suppose that 
F is non-tangentially bounded at every point of a set £■ C M. Then 

F{uj + iO) = lim F{uj + ir;) (6.3) 
exists for almost every lu (z E. 

Note: In fact, the so-called non-tangential limits lim^^j^ with C re- 

stricted to re(cij), exist at almost every uj, but we will not use this fact. 
Our principle theorem on losses is the following 

Theorem 6.2 Suppose the mass operator m is strictly positive, m > 6 > 
0, and the Laplace transform d (C) of the friction function satisfies the 
following strengthened form of the power-dissipation condition (|2.42|) : there 
are measurable functions 6 : M. ^ (0,7r), 7 : R — (0,oo), with 7(-) S 
L;^Q^(M) such that for almost every w G R 

Re{v,a(C)v) > for all C £ re(,^)(u;) and w € (6.4) 

Then for any compactly supported generalized force / e Ll{M.;Ho), the 
solution Vf (t) to the evolution equation Hl.l|l vanishes in the large t limit: 

hm \\vj{t)\\ = 0. (6.5) 

Remark 6.3 The condition H6.4|l holds in particular if Red(C) is strictly 

positive, Rea(C) > SIhq for all ImC > 0. For example, in the Lorentz 

2 

medium we have a (C) = 7 + i^, so Re a {() > jIho- 

Proof. Note that the admittance 21 (C) = i (C"^ — A + ia (C))~^ obeys 

Re2t(C) =2t(C) [ImCm + Rea(C)]2t(C)^ 

> ^ 21 (C) 21 (C)^ , for C e r,(^) {u;). (6.6) 



7(w) 
We conclude that 



1121(011 > |lRc2l(C)|l >—- |12l(C)f, (6.7) 
7(cj) 
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and thus 

1121(011 < 7H, force re(^)(a;). (6.8) 

Using Theorem 13.131 we find a conservative extension consisting of a 
self-adjoint operator ^ on a Hilbert space TC and an isometric imbedding 

: Ho^n such that 21 (C) = im-^/^T (C - ^1)"^ T^m,-'^/'^. Furthermore 
the solution Vf to with generalized force / is 

Vf{t)= I m-^^^Te'^''-'^-^T^m-^/^f{s)ds. (6.9) 

The theorem now follows from Theorem 16.41 below via dominated conver- 
gence. ■ 

Theorem 6.4 Let ^ be a self adjoint operator on a separable Hilbert 
space Ti., and let T : 7i ^ Hq be a bounded map from Ti. into a separable 
Hilbert space Hq. Suppose there are measurable functions 7 : R — > (0, 00) 
and 6* : K ^ (0, tt) with 7(-) e Lj^^iM.) such that for almost every w € IK 

sup ||lmr(C-^)"^r^|| < 7(w). (6.10) 

Then 

lim Te""-^T^w = (6.11) 

t — ►iboo 

for every v G Hq. 

Remark 6.5 (1) As will be clear from the proof, if 7(-) G (instead of 
-^ioc)' then eq. H6.11|l holds in norm. (2) Eq. (|6.10() implies that the 
spectral measure associated to A and any ip G ranT^^ is purely absolutely 
continuous, with density 

- lim (?/>,Im(tj + iry-yl)"V) ■ (6.12) 

e— ++0 

Proof. We begin by showing that for each bounded open interval / C M 

||Te-'*-^£;^(/)r^|| ^ : as : t ±00, (6.13) 

where Ea{-) is the spectral resolution of A. 
Let us define 

F{C) = ImT^(C - A)-^T for ImC ^ 0. (6.14) 

By the spectral theorem and the observation that the spectral measures of 
A, d\\E_A{X)i/j\\'^ , are purely absolutely continuous, we have 

Te-''^EAiI)T^v = --lim [ da; e-'''^ F {lj + hj) v , (6.15) 

TT riia J J 
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for each finite t and every v G Hq 

Since Hq is separable, it has a countable basis (pj, j — f,....,. Using 
Theorem 16. II we conclude that there is /' C / with m{I\I') — such that 

lim(0„F(w + i77)0,) (6.16) 

exists for every oj G I' and every pair of basis vectors (pi , (pj ■ Let S denote 
the subspace of linear combinations of finitely many basis vectors. Given 

oo oo 

Vlim \{<f>,,F{Lj + ir])i:}\'^ < liminf V | F(w + ir/)V') |^ 

^ — ^ rj^O 77^0 — ^ 



i=l 



= liminf ||F(w + i?7)V'|r < 7('^)^||V'II^- (6.17) 
We conclude that, for every uj G I' , the limit 

wk-limF(cj + ifj)^ = V \im{4:,,F{uj + i'n)^p)(p^ (6.18) 

exists, and the map 



is bounded from S into Hq. Since S is dense, this map may be extended to a 
unique bounded linear map i^(a; + iO) : Hq Hq with \\F{iLj + iO)|| < j{u). 
It is elementary to see that, 

(V-i, F(w + 10)^^2) = lim {^pi,F{uj + bj)ip2) for any ^2 e iJo, (6.20) 
77— >o 

i.e. 

F{lj + iO) = wklimF(w + i?^) for uj e /'. (6.21) 
By dominated convergence we find from (|6.15|) 



I" J I 

for every v G Ha- Since + iO)|| < 7(0-1) G L^{I), we have 



(6.22) 



lim ||re-'*-^£;^(/)r1| = (6.23) 

by the Riemann-Lebesgue lemma - the extension of this result to operator 
valued functions is elementary. 
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To complete the proof, we note that given e > and G Hq we can 
find a finite interval 4 C M such that ||(1 - £;^(/e))TV|| < £• Thus 

limsup ||Te""-^rV|| < limsup ||re""-^£;^(/e)rV|| + e = £• (6.24) 

t — ^±oo t — >±oo 

Since e is arbitrary, we see that H6.11|l holds. ■ 

7 Operator versions of classical spectral the- 
orems 

In this section we discuss the proofs of operator versions of Bochner's The- 
orem 13.11 - Theorems 13.21 and 13.61 above - and operator versions of the 
Herglotz-Nevanlinna Theorems IT^ 13 . 81 - Theorcm l3. 13l and l3. 161 For prop- 
erties of operator-valued functions holomorphic in a half-plane and their 
boundary values see [T^ . 

7.1 Bochner's Theorem 

The two operator valued generalizations of Bochner's theorem stated above 
- Theorems 13 . 21 and 13 . 61 - are combined in the following statement. 

Theorem 7.1 The friction function Ue {t) = 2aoo5 (t) + ae{t), — oo < t < 
oo, with ae{t) a strongly continuous B {Hq) valued function, and aoo a 
bounded non-negative operator is representable as 

ae (t) = Dhmre-""^ {T<i>%y , $1 = + ^H.) \ (7.1) 

with ill a self-adjoint operator on Hi and F : D (fli) Hq a fii-bounded 
linear map, if and only if a (t) satisfies the dissipation condition (|1.19() 
for every continuous Hq valued function v{t) with compact support. The 
operator F is bounded if and only if aoo = 0, in which case 

ae(t) = re~'*"iF^ (7.2) 

If the space Hi is minimal - in the sense that 

((F/(r!i))%: /eC,(M), veHo) (7.3) 

is dense in Hi - then the triplet {Hi, r2i,F} is determined uniquely up to 
an isomorphism. 
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Proof. Let us start by defining the Hilbert space Hi, which we take 
to be an extension of the Banach space V — (R, Hq) n (R, Hq) of 
measurable Hq valued functions with 



1/2 

2 



Um'dt] + |10(OMi<oo. (7.4) 

i / J — oo 

Let (•, denote the inner product on Hq and define on ^ a quadratic 
form 



oo poo 



(</>, 1^) := / / (0 (t) ,aAt-T)i; (r))^^ dtdr (7.5) 

OO J — OO 

(<^(t),aooV'(i))H„ (7.6) 
(0 {t) ,ae{t-T)il) (r))^^^ didr, 



— oo 

oo /"OO 



oo ^ — oo 



which is positive semi-definite by virtue of the power dissipation condition. 
However, there may be null vectors in V , that is vectors (j) with (0, (p) ~ 0. 
Let N denote the set of null vectors and define Hi = V/N , where " " " 
denotes closure in the norm inherited from the inner product. Then Hi is 
a Hilbert space, whose inner product we also denote by (•,•). 

Any <j> G V defines a unique element [(p] G V/N C Hi with = 
{(f>, (t>). In particular, [p] is zero if and only if (0, (p) = 0. Furthermore, the 
map p [pi] is bounded from V into Hi since 



/oo / poo 

\\P{t)f dt + sup||ae(r)|| / 
-oo •'■ \J — oo 



2 

0(011 dt ' 



<ml- (7.7) 



Thus a convergent sequence pj in V gives rise to a convergent sequence 
[(pj] in Hi and 



lim [0j 



lim pj 

j 



We define the operator fii to be the self adjoint generator of the one 
parameter unitary group of time translations. Specifically, we note that 
the transformations Tg of V into itself given by 

T,0(i) =0(t-s), (7.9) 

form a group which preserves the pre- inner product H7.5() . Therefore, s i— > 
Ts extends to a one parameter unitary group s i— s- Us on Hi , which is in fact 
strongly continuous (as follows from strong continuity of on V). The 
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Stone-von Neumann theorem implies there is a unique self adjoint operator 
fii on Hi with Us = e'^^^i . 

Clearly there is a connection between fii and the operation of differ- 
entiation on y. To understand this, note that given € we have 
[(j)] e ^(r^i) if and only if 



[0] ^ hm ^MtM = _i lin, 
s^O is s— >o 



(7.10) 



exists. Thus a sufficient condition for [(j)] to be in V (fii) is for s ^ {Tgcf) — (p) 
to converge in V, which holds if and only if ^ S D (9*), in which case 

lim^^ = -9.0. (7.11) 

s S 

Therefore we have 

Qi [(j)]^i[dt(j)] for cj) e V (dt) . (7.12) 

Note, however, that [(/)] e V (fli) does not necessarily imply e I? (9t). 

If we formally define to be the map T^v = [v6 (•)], ignoring for the 
moment that vS {■) ^ V, we may calculate that 



ae{-t)(j){t) dt, (7.13) 

J — OO 

and therefore 

POO 

rg-itfiipt^; = / ae(-s)w(5(s + <) = ae(t)w. (7.14) 



If aoo = then in fact H7.13|l defines a bounded operator, and the above 
calculation may be justified. In that case, we could have started with a 
space V including point measures so that vS (t) for v G Hq would be in V 
and we would have T^v — [vS (•)]. 

However, to consider also aoo 0, we work indirectly by defining the 
bounded map 5^ ; Ho Hi, 

S^v := [vG] ; G (t) = ie* J ^ ^ ° ' (7.15) 
^ ^ ' ^ ^ [0 i > . ^ ' 

Note that formally, {Qi ~ Uh,)S^v = [v6 (•)], since {idt ~i}G{t) ^6 (t). 
To proceed rigorously, let us compute S := ("S*^)^. Given (j) e V 



oo />oo 



{v,Sm = {S^v,[^])^ / G{t){v,ae(t-T)cl){T)) dtdT, (7.16) 



oo -/ — oo 
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so 



^-i J e*ae {t - t) d?j (r) dr 

-i / e^-*ae i-t) dt ] 4> (r) dr. (7.17) 



If we define V : V [VLi) ^ Hq hy 

r := 5 + i/ffj , (7.18) 
then we recover (|7.13|) for (f) £ T> (dt), 

oo / />oo 

T-t 



— OO \>y r 



e^-'fle (-t) dtj {dr + 1} (r) dr 

POO 

fle (-r) (t) dr = 2aoo9i (0) + / (-r) (t) dr, (7.19) 



since 4> G 'D {dt) implies (/) is continuous (so (0) is unambiguous). 
One may easily verify that 

= e-l*lao, + -/ e-l*+^lae(s) ds, (7.20) 

and therefore 

^ + l) ^e"'"'^ 5^ = 2aoo<5 (i) + (t) = (t) . (7.21) 



dt2 ^ 

Thus, with $7? defined by H7.1|l and using (|7.18|1 and H7.20|l we get 



re--(r<.|)t.(4 + i)r^e--(r. 



2Nt 

r) 

+ 1^ Se-''"'^lS^ ^ ae (t) as i? ^ oo, (7.22) 

where —> denotes limit in the sense of distributions. 

The uniqueness up to isomorphism can be understood as follows. Let 
{i?i,f2i,r} and {H'i,il'i,r'} be distinct representations, and suppose that 

S = ({Tf (f^i))^ v. feC,{R),ve H„) (7.23) 
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is dense in Hi . We denote by / the Fourier transform of / G Cc (K) , so 
that 

/oo 
f{t)e'''''dt. (7.24) 
-OC 

Then, given f.g^Cc (M) and v,w ^ Hq, we see that 

I H 

: lim (w,Tg{ni)f{ni)^ {T^l)\ 

R-^oc \ 

/oc nOC 
/ g (t) f* (r) (w, re'(^-*)"i (r^^)^ v) dtdr 

OO nOC \ ' --^^ ) 

/ 5 (0 /* ('t) ("^i (i - t) w) didr 

-OO •/ — OO 

= ((r'g(f^;))^z«,(r'/(f^'i))^)^,. 

Thus, defining 

r(r/(f^i))^t; = (r'/(f^'i))^ (7.26) 

and extending T to 5 by hnearity, we produce a well defined isometry 
T : S ^ H[. The closure of this map, also denoted T, is an isometric 
imbedding T : Hi '-^ H[. It is now easy to verify that 

/ {ni) = rt/ (r!'i) T, for / e (K) ; r' = rrt, (7.27) 

and thus the representation {i?i,f2i,r} is isomorphic to the restriction of 
{H'l.n'iX} to the closure of 

((r7 {n'i) f v. / e (M) , « e i/o) . (7.28) 



7.2 Herglotz-Nevanlinna Theorems 

In Section 2, we have already presented a more or less complete proof 
of Theorem 13.131 based on the classical Nevanlinna Theorem 13.81 and the 
Naimark representation ITT^ of generalized spectral families. Since Theo- 
rem |^^] is proved in a very similar fashion, we present only a somewhat 
streamlined proof here. 

Proof of Theorem 13.161 By the Herglotz-Nevanlinna theorem, for 
each V G Hq there is a finite measure dNy^y and a real number S,v,v such 
that 

/oo -1 I > 
—^dN,,, (a) . (7.29) 
-oo cr-C 
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The 9 term drops out of the representation because ( {() ^ as ( ^ 
oo. In fact, ^y^y = Re (w, G (i) w). As in the derivation of Theorem 
13.131 we define "off-diagonal" measures dNy^w for each pair v,w G Hq by 
polarization so that 

/c=o -1 I 
—^dNy^y, (a) . (7.30) 



As above there is a "generalized spectral family" K {a) satisfying the 
hypothesis of the Naimark Theorem 13.121 such that 



{v,K{a)w) - / dNy^y,{cj). (7.31) 

J ( — oOjCr] 

We denote by Hi and E (cr) the Hilbert space and resolution of the iden- 
tity guaranteed by the Naimark Theorem, letting T denote the associated 
mapping T : Hi ^ Hg. Thus, we have 

G(C) =ReG(i)+ ii^Tdi;(CT)rt, (7.32) 



or 



G(C)=ReG(i)+T^f^±^rt, r adE{a). (7.33) 

"1 - MHi J-oo 

From this we may easily compute the following formula 



Therefore, given v £ Hq, 



(7.34) 



G{Ov= Um r^^ — ^ ( T-^^j—^ ) , r = t^Ih, + nl (7.35) 



^ fr^' 

fl^oo Qi - CIh, \ R^Ih, + 
because we have 



lim {ReG(ii?)}?; = (7.36) 



by assumption, and 



lini T-iM^-J^l±^Tt. = (7.37) 
for each ( in the upper half plane, since ^^^^^^^^^ — > strongly. ■ 
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A Appendix: Stieltjes Inversion Formula and 
Naimark's Theorem 



In Section 3 - and also in the examples of Section 4 - we discussed the 
construction of conservative extensions. Reference was made there to two 
classical constructions - the Stieltjes Inversion formula and Naimark's con- 
struction for generalized spectral measures - which quite generally provide 
an explicit description of the Hilbert space Hi in the operator versions of 
the Herglotz-Nevanlinna theorems. For completeness we include a discus- 
sion of those results here. 



A.l Stieltjes Inversion Formula 

The Nevanlinna Theorem 13. 81 - in particular the relation H3.15|l - suggests 
the introduction of the so-called Cauchy transform defined for complex- 
valued measures of finite variation on M: 

r^.:„,C#o. (A.1) 

The Nevanlinna theorem states that the set of functions which are Cauchy 
transforms of non-negative finite measures is exactly the class of analytic 
maps of the upper half plane into itself which decay as O (1/ Im as Im C, 
oo. There does not seem to be such a simple description of the set of 
Cauchy transforms of complex measures. Nonetheless, a complex measure 
is uniquely determined by its Cauchy transform, JOl (Section 32.1, Lemma 
4)], 

Proposition A.l The Cauchy transform (|A.ip is one-to-one, i.e. a com- 
plex measure of finite variation is uniquely determined by its Cauchy trans- 
formation. Furthermore, if dN (a) is a real (signed) measure of finite vari- 
ation it can be recovered from its Cauchy transform TV (C) restricted to 
{ImC > 0} by Stieltjes' formula : 

AT (ai + 0) + TV (fJi - 0) _ A^ (gp + 0) + A^ ((Tq - 0) ^ 

2 2 ^ ' ' 



1 

= lim — / Im N (a m) da. 



If n (ct) da is the absolutely continuous component of the measure dN (a), 
in particular if dN (a) = n (a) da, we also have 

n (a) = lim — ImTV {a + irj) for Lebesgue almost every a. (A. 3) 
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Remark A. 2 Another manifestation of ljA.2|l . is the weak convergence 



wkhm - Im N (cr + iij) da = dN (cr) , (A. 4) 



that is, 

/(a) dN{<j) 



1 f°° 

= hm - / / (cr) lmN{a + ir;) da for all f e Cc (M) . (A.5) 

Thus a scalar Herglotz function 5 {() which is 0(1/ Im () as Im C — > 00 
may be represented by the formula 

1 1 

5 (C) = lim - / lmgia + irj) da. (A.6) 

In particular, if Img (ct + iO) ~ lim,,^o Img (cr + irj) exists for almost every 
a and 

/oo 
|Im g (cr + IT]) - Im g (cr + iO) | dcr = 0, (A. 7) 

-00 



then ^ 

g{0 = - I Im .g {a + iO) da. (A.8) 

In general, however there may be a singular component to the measure dN 
in the representation H3.15|l . 

For an operator valued Herglotz function G (C), with \\G (C)|| = O (1/ ImC), 
we have therefore 

G (C) = lim - f Im G (ct + irj) da, (A.9) 

'f^o TT J^^ a-Q 

with the integral understood in the weak sense, i.e. 
{v,G{Ow) 

1 [°° 1 

= lim — / (w, Im G (cr + 177) w) da, for v,w & Hq. (A. 10) 

77^0 TT J^^ a-Q 

Thus the generalized spectral family K [a) associated to G (C) can be ex- 
pressed through the formula 

d (u, K (cr) w) — wklim — (w, Im G (cr + iry) w) da. (A. 11) 
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A. 2 Naimark's Theorem 



The Naimark construction for K (cr), which leads to Theorem 13. 121 is most 
easily understood by realizing Hi as the Hilbert space (dK), where the 
latter space needs to be appropriately defined. Formally (dK) should 
consist of all Hq valued functions 4" such that (\E' (cr) , dK (a) ^' (tr)) is 
finite, modulo null functions (for which the integral is 0). It is not always 
clear how to make sense of this integral however, and one must turn to a 
more abstract definition of (K). However, if dK (a) = m (a) da, with 
m (cr) bounded for almost every cr, we can avoid the abstract construction 
by defining Hi ~ L? {dK) to be the space of Hq valued functions such 
that 

/oo 
(«'(cr),m(cr)^'(cr)) da<oo, (A.12) 
-oo 

modulo null functions for which ||^||^ — 0. Note that m (cr) is necessarily 
a positive operator since 



1 



cr+e 



(w, m (ct) w) = lim — / d{v,K{(j)v). (A.13) 

For more general measures dK there are two options. We could define 
the norm appearing in (|A.12|1 for Hq valued simple functions - functions 
taking a finite number of values - and let Hi be the closure of the space of 
simple functions under this norm. Alternatively, we can express dK (cr) as 
the weak limit (|A.11|I and define Hi — (dK) as the space of functions 
\I> such that 

ll^-f -lim-/ (^((j),lmG(a + iT])^(a)) da (A.14) 

exists and is finite, modulo null functions (as always.) 

However it may be defined, once the space Hi — I? (dK) has been 
constructed, there is a natural spectral measure E (cr) given by 

[E{am{-) = ^ if.;. ' (^-i^) 

that is E (cr) corresponds to multiplication by the characteristic function of 
(— oo, cr]. The associated self adjoint operator Vli = adE (cr) is simply 
multiplication by the independent variable: 

fli^ (a) ^ a^ {a) . (A.16) 

Finally there is a natural map F^ : Hq Hi which takes an element 
ip G Hq to the constant function with value ip: 

[T^] (cr) = V for every cr G M. (A.17) 
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It is easy to verify that 

K{a) =TE{a)T^. (A. 18) 
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